This code is Script 4a for Kitchel et al. Biotic homogenization, the exception and not the rule for marine fish communities manuscript.

COPY FROM THIS SCRIPT

(We expect this to be the main analyses in the text)

SESSION INFO TO DO

library(data.table)
library(vegan)
library(sf)
library(concaveman) #polygon around points
library(betapart) #allows us to partition beta diversity
library(geosphere)
library(ggpubr) #stat_regline_equation
library(nlme)
library(mgcv) #to make gam
library(cowplot)
library(lme4)

#Pull Dissimilarity Means
distances_dissimilarities_allyears.r <- readRDS(here::here("output", "dissimilarities", "distances_dissimilarities_allyears.r.rds"))

#make survey and survey unit factors
distances_dissimilarities_allyears.r[,survey:=factor(survey)][,survey_unit:=factor(survey_unit)]

#adjust years
distances_dissimilarities_allyears.r[,year_adj := year-min(year)+1]

#add new variable for year in sequence per region
distances_dissimilarities_allyears.r[,first_year := min(year),.(survey_unit)]
distances_dissimilarities_allyears.r[,last_year := max(year),.(survey_unit)]

#distances_dissimilarities_allyears.r[,year_in_seq := year-first_year+1]

distances_dissimilarities_allyears.r[,years_sampled := last_year-first_year+1]

###Palette for Plotting Palette for plotting all 42 survey units

survey_unit.list <- levels(factor(distances_dissimilarities_allyears.r[,survey_unit]))

palette_42 <- c(
  "#5A5156", #AI
  "#DF00DB", #BITS-1
  "#DB8EDA", #BITS-4
  "#F6222E", #CHL
  "#F8A19F", #DFO-NF
  "#16FF32", #DFO-QCS
  "#325A9B", #EBS
  "#3283FE", #EVHOE
  "#FEAF16", #FR-CGFS
  "#fccb6d", #GMEX-Fall
  "#1C8356", #GMEX-Summer
  "#C4451C", #GOA
  "#85660D", #GRL-DE
  "#B0009F", #GSL-N
  "#BF79B8", #GSL-S
  "#1CBE4F", #ICE-GFS
  "#782AB6", #IE-IGFS
  "#90AD1C", #MEDITS
  "#6B003A", #NAM
  "#A75B00", #NEUS-Fall
  "#E3B072", #NEUS-Spring
  "#02E8B6", #NIGFS-1
  "#97E7D5", #NIGFS-4
  "#B00068", #Nor-BTS-3
  "#00B9E3", #NS-IBTS-1
  "#95E2F4", #NS-IBTS-3
  "#B3CE73", #NZ-CHAT
  "#689500", #NZ-ECSI
  "#364d02",#NZ-SUBA
  "#AAF400", #NZ-WCSI
  "#AA0DFE", #PT-IBTS
  "#7f9eb8", #ROCKALL
  "#FA0087", #S-GEORG
  "#DEA0FD", #SCS-Summer
  "#FCEF88", #SEUS-fall
  "#A59405", #SEUS-spring
  "#FCE100", #SEUS-summer
  "#544563", #SWC-IBTS-1
  "#a37fc7", #SWC-IBTS-4
  "#C075A6", #WCANN
  "#BDCDFF", #ZAF-ATL
  "#003EFF"  #ZAF-IND
)

color_link <- data.table(survey_unit = survey_unit.list,hex = palette_42)

Add names for plotting


name_helper <- data.table(Survey_Name_Season = c("Aleutian Islands",
                                    "Baltic Sea Q1",
                                    "Baltic Sea Q4",
                                    "Chile",
                                    "Newfoundland",
                                    "Queen Charlotte Sound",
                                    "Eastern Bering Sea",
                                    "Bay of Biscay",
                                    "English Channel",
                                    "Gulf of Mexico Summer",
                                    "Gulf of Alaska",
                                    "Greenland",
                                    "N Gulf of St. Lawrence",
                                    "S Gulf of St. Lawrence",
                                    "Iceland",
                                    "Irish Sea",
                                    "Mediterranean",
                                    "Namibia",
                                    "NE US Fall",
                                    "NE US Spring",
                                    "N Ireland Q1",
                                    "N Ireland Q4",
                                    "Barents Sea Norway Q3",
                                    "N Sea Q1",
                                    "N Sea Q3",
                                    "Chatham Rise NZ",
                                    "E Coast S Island NZ",
                                    "W Coast S Island NZ",
                                    "Portugal",
                                    "S Georgia",
                                  "Scotian Shelf Summer",
                                  "SE US Fall",
                                  "SE US Spring",
                                  "SE US Summer",
                                  "W Coast US",
                                  "Atlantic Ocean ZA",
                                  "Indian Ocean ZA",
                                   "Rockall Plateau",
                                  "Scotland Shelf Sea Q1",
                                  "Scotland Shelf Sea Q4",
                                  "Falkland Islands",
                                  "Gulf of Mexico Fall",
                                  "Barents Sea Norway Q1",
                                  "Sub-Arctic NZ",
                                  "Scotian Shelf Spring"),
                          survey_unit = c(
                                  "AI",        
                                  "BITS-1",    
                                  "BITS-4",    
                                  "CHL",       
                                  "DFO-NF",    
                                  "DFO-QCS",   
                                  "EBS",       
                                  "EVHOE",     
                                  "FR-CGFS",   
                                  "GMEX-Summer",
                                  "GOA",       
                                  "GRL-DE",    
                                  "GSL-N",     
                                  "GSL-S",     
                                  "ICE-GFS",   
                                  "IE-IGFS",   
                                  "MEDITS",    
                                  "NAM",       
                                  "NEUS-Fall", 
                                  "NEUS-Spring",
                                  "NIGFS-1",   
                                  "NIGFS-4",   
                                  "Nor-BTS-3", 
                                  "NS-IBTS-1", 
                                  "NS-IBTS-3", 
                                  "NZ-CHAT",   
                                  "NZ-ECSI",   
                                  "NZ-WCSI",   
                                  "PT-IBTS",   
                                  "S-GEORG",   
                                  "SCS-SUMMER",
                                  "SEUS-fall", 
                                  "SEUS-spring",
                                  "SEUS-summer",
                                  "WCANN",     
                                  "ZAF-ATL",   
                                  "ZAF-IND",
                                  "ROCKALL",
                                  "SWC-IBTS-1",
                                  "SWC-IBTS-4",
                                  "FALK",
                                  "GMEX-Fall",
                                  "Nor-BTS-1",
                                  "NZ-SUBA",
                                  "SCS-SPRING"
                          ))


color_link <- name_helper[color_link, on = "survey_unit"]

##Make GAMs

Bray Curtis

bray_curtis_total_gam <- gam(bray_curtis_dissimilarity_total_mean ~ year + s(survey_unit, year, bs = "fs", m = 1),#random smooth
                            data = distances_dissimilarities_allyears.r)

##Make LMERS

Bray These all converged

#running with lme instead of lmer gave same results, but allowed for calculation of p-value
bray_curtis_total_lme <- lme(bray_curtis_dissimilarity_total_mean ~ year_adj, random = (~1 + year_adj|survey_unit),data = distances_dissimilarities_allyears.r)

#but also run with lmer for confint
#total
bray_curtis_total_lmer <- lmer(bray_curtis_dissimilarity_total_mean ~ year_adj + (1 + year_adj|survey_unit),data = distances_dissimilarities_allyears.r)
Warning: Model failed to converge with max|grad| = 0.137452 (tol = 0.002, component 1)
#balanced (for comparison plot)
bray_curtis_balanced_lmer <- lmer(bray_curtis_dissimilarity_balanced_mean ~ year_adj + (1 + year_adj|survey_unit),data = distances_dissimilarities_allyears.r)
#gradient (for comparison plot)
bray_curtis_gradient_lmer <- lmer(bray_curtis_dissimilarity_gradient_mean ~ year_adj + (1 + year_adj|survey_unit),data = distances_dissimilarities_allyears.r)
Warning: Model failed to converge with max|grad| = 0.012616 (tol = 0.002, component 1)
#Jaccard total (for comparison plot)
jaccard_total_lmer <- lmer(jaccard_dissimilarity_total_mean ~ year_adj + (1 + year_adj|survey_unit),data = distances_dissimilarities_allyears.r)
Warning: Model failed to converge with max|grad| = 0.00568546 (tol = 0.002, component 1)
summary(bray_curtis_total_lme)
Linear mixed-effects model fit by REML
  Data: distances_dissimilarities_allyears.r 

Random effects:
 Formula: ~1 + year_adj | survey_unit
 Structure: General positive-definite, Log-Cholesky parametrization
            StdDev      Corr  
(Intercept) 0.108413092 (Intr)
year_adj    0.001972694 -0.763
Residual    0.023981944       

Fixed effects:  bray_curtis_dissimilarity_total_mean ~ year_adj 
 Correlation: 
         (Intr)
year_adj -0.791

Standardized Within-Group Residuals:
       Min         Q1        Med         Q3        Max 
-4.8791736 -0.4684941  0.0317567  0.5158154  4.3793155 

Number of Observations: 897
Number of Groups: 42 
anova(bray_curtis_total_lme)

bray_curtis_total_coefs <- data.table(coef(bray_curtis_total_lme))
bray_curtis_total_coefs[,survey_unit := rownames(coef(bray_curtis_total_lme))][,Year := round(year_adj,5)][,Intercept := round(`(Intercept)`,2)]
#View(bray_curtis_total_coefs)

bray_curtis_total_coefs <- bray_curtis_total_coefs[color_link, on = "survey_unit"]

bray_curtis_total_coefs.exp <- bray_curtis_total_coefs[,.(Survey_Name_Season, Intercept, Year)]

#export this table
fwrite(bray_curtis_total_coefs.exp, file = here::here("output","bray_curtis_total_coefs.exp.csv"))

Get LMER model as predictions


# need to sort out year in seq versus overall year models
#new data for lmer
lmer_year <- seq(min(distances_dissimilarities_allyears.r[,year]), max(distances_dissimilarities_allyears.r[,year]), by = 0.1)

lmer_year_adj <- seq(min(distances_dissimilarities_allyears.r[,year_adj]), max(distances_dissimilarities_allyears.r[,year_adj]), by = 0.1)

#predict average lmer
lmer_bray_total_predictions <- data.table(year = lmer_year, year_adj = lmer_year_adj)

#confidence intervals
bray_curtis_total_lmer_confint <- confint(bray_curtis_total_lmer)
Computing profile confidence intervals ...
#populate data table of lmer predictions
lmer_bray_total_predictions[,bray_curtis_lmer_preds := fixef(bray_curtis_total_lmer)[[1]] + fixef(bray_curtis_total_lmer)[[2]] * year_adj][,bray_curtis_lmer_preds_lowerCI := bray_curtis_total_lmer_confint[5] + bray_curtis_total_lmer_confint[6] * year_adj][,bray_curtis_lmer_preds_upperCI := bray_curtis_total_lmer_confint[11] + bray_curtis_total_lmer_confint[12] * year_adj]

# need to sort out year in seq versus overall year models
#new data for lmer


#predict average lmer
lmer_bray_balanced_predictions <- data.table(year = lmer_year, year_adj = lmer_year_adj)

#confidence intervals
bray_curtis_balanced_lmer_confint <- confint(bray_curtis_balanced_lmer)
Computing profile confidence intervals ...
#populate data table of lmer predictions
lmer_bray_balanced_predictions[,bray_curtis_lmer_preds := fixef(bray_curtis_balanced_lmer)[[1]] + fixef(bray_curtis_balanced_lmer)[[2]] * year_adj][,bray_curtis_lmer_preds_lowerCI := bray_curtis_balanced_lmer_confint[5] + bray_curtis_balanced_lmer_confint[6] * year_adj][,bray_curtis_lmer_preds_upperCI := bray_curtis_balanced_lmer_confint[11] + bray_curtis_balanced_lmer_confint[12] * year_adj]

# need to sort out year in seq versus overall year models


#predict average lmer
lmer_bray_gradient_predictions <- data.table(year = lmer_year, year_adj = lmer_year_adj)

#confidence intervals
bray_curtis_gradient_lmer_confint <- confint(bray_curtis_gradient_lmer)
Computing profile confidence intervals ...
#populate data table of lmer predictions
lmer_bray_gradient_predictions[,bray_curtis_lmer_preds := fixef(bray_curtis_gradient_lmer)[[1]] + fixef(bray_curtis_gradient_lmer)[[2]] * year_adj][,bray_curtis_lmer_preds_lowerCI := bray_curtis_gradient_lmer_confint[5] + bray_curtis_gradient_lmer_confint[6] * year_adj][,bray_curtis_lmer_preds_upperCI := bray_curtis_gradient_lmer_confint[11] + bray_curtis_gradient_lmer_confint[12] * year_adj]

# need to sort out year in seq versus overall year models


#predict average lmer
lmer_jaccard_total_predictions <- data.table(year = lmer_year, year_adj = lmer_year_adj)

#confidence intervals
jaccard_total_lmer_confint <- confint(jaccard_total_lmer)
Computing profile confidence intervals ...
#populate data table of lmer predictions
lmer_jaccard_total_predictions[,jaccard_lmer_preds := fixef(jaccard_total_lmer)[[1]] + fixef(jaccard_total_lmer)[[2]] * year_adj][,jaccard_lmer_preds_lowerCI := jaccard_total_lmer_confint[5] + jaccard_total_lmer_confint[6] * year_adj][,jaccard_lmer_preds_upperCI := jaccard_total_lmer_confint[11] + jaccard_total_lmer_confint[12] * year_adj]

###Move forward with Bray Curtis total (others to supplement, although will include in plot below)

Coefficients for LMER by survey_unit

#unique survey unit year combos
survey_unit_sampling_years <- unique(distances_dissimilarities_allyears.r[,.(survey_unit, year_adj, year, years_sampled)])

# see group coefficients
model_coefs_reduced <- data.table(transform(as.data.frame(ranef(bray_curtis_total_lmer)), lwr = condval - 1.96*condsd, upr = condval + 1.96*condsd))
#https://stackoverflow.com/questions/69805532/extract-the-confidence-intervals-of-lmer-random-effects-plotted-with-dotplotra


#ONLY SLOPES
model_coefs_reduced <- model_coefs_reduced[term == "year_adj",]

model_coefs_reduced[,survey_unit := grp][,year_adj := condval]

#merge with duration of survey
model_coefs_reduced_length <- model_coefs_reduced[survey_unit_sampling_years, on = "survey_unit"]


model_coefs_reduced_length[,years_sampled := as.numeric(years_sampled)][,Directional_Change := ifelse(year_adj > 0, "Differentiation","Homogenization")]

#does it cross zero?
model_coefs_reduced_length[,significant := ifelse(lwr >0 & upr>0,T,ifelse(lwr<0 & upr<0,T,F))]

#significant directional change
model_coefs_reduced_length[,Directional_Change_sig := ifelse(year_adj > 0 & significant == T, "Differentiation",ifelse(year_adj < 0 & significant == T, "Homogenization", "No trend in\ndissimilarity"))]


#min max distances_dissimilarities
min_bray_reduced <- min(distances_dissimilarities_allyears.r[,bray_curtis_dissimilarity_total_mean], na.rm = T)
max_bray_reduced <- max(distances_dissimilarities_allyears.r[,bray_curtis_dissimilarity_total_mean], na.rm = T)

model_coefs_reduced_length <- model_coefs_reduced_length[color_link, on = "survey_unit"]

#delete any NAs
model_coefs_reduced_length <- na.omit(model_coefs_reduced_length, cols = "significant")

#order table by coefficient
setorder(model_coefs_reduced_length, year_adj)

BC_total_model_coefs_reduced_length.unique <- unique(model_coefs_reduced_length[,.(condval,condsd, lwr, upr, survey_unit, year_adj, years_sampled, Directional_Change, hex, Survey_Name_Season, significant, Directional_Change_sig)]) 

#extract color hexes
#year adj coef order
color_year_adj_order <- BC_total_model_coefs_reduced_length.unique[,hex]

#alphabetical order
BC_total_model_coefs_reduced_length.unique.alpha <- setorder(BC_total_model_coefs_reduced_length.unique, Survey_Name_Season)

BC_total_model_coefs_reduced_length.unique.alpha[,trend_color := ifelse(Directional_Change_sig == "Homogenization", "#e7ac5b", ifelse(Directional_Change_sig == "Differentiation","#91c874","#cbbfde"))]

color_alpha_order <- BC_total_model_coefs_reduced_length.unique.alpha[,hex]
color_alpha_order_bytrend <- BC_total_model_coefs_reduced_length.unique.alpha[, trend_color]

saveRDS(BC_total_model_coefs_reduced_length.unique, here::here("output","region_stats","BC_total_model_coefs_reduced_length.unique.Rds"))

Total versus balanced BC plot

ggplot(distances_dissimilarities_allyears.r) +
  geom_point(aes(bray_curtis_dissimilarity_total_mean, bray_curtis_dissimilarity_balanced_mean)) +
  geom_abline(slope = 1, intercpet = 0) +
  geom_smooth(aes(bray_curtis_dissimilarity_total_mean, bray_curtis_dissimilarity_balanced_mean)) +
  theme_classic() +
  labs(x = "Total BC dissimilarity",  y = "Balanced changes in abundance/biomass")
Warning: Ignoring unknown parameters: `intercpet`

ggplot(distances_dissimilarities_allyears.r) +
  geom_point(aes(bray_curtis_dissimilarity_total_mean, bray_curtis_dissimilarity_gradient_mean)) +
  geom_abline(slope = 1, intercpet = 0) +
  geom_smooth(aes(bray_curtis_dissimilarity_total_mean, bray_curtis_dissimilarity_gradient_mean)) +
  theme_classic() +
  labs(x = "Total BC dissimilarity",  y = "Abundance/biomass gradient")
Warning: Ignoring unknown parameters: `intercpet`
#all on one plot
#reduced
distances_dissimilarities_allyears.r.r <- distances_dissimilarities_allyears.r[,.(year, bray_curtis_dissimilarity_total_mean, bray_curtis_dissimilarity_gradient_mean, bray_curtis_dissimilarity_balanced_mean, jaccard_dissimilarity_total_mean)]

#wide to long
dissim_long <- melt(distances_dissimilarities_allyears.r.r, id.vars = c("year"), variable.name = "Dissimilarity metric")

dissim_long[,`Dissimilarity metric` := factor(`Dissimilarity metric`, levels = c( "bray_curtis_dissimilarity_total_mean","bray_curtis_dissimilarity_balanced_mean","bray_curtis_dissimilarity_gradient_mean", "jaccard_dissimilarity_total_mean"))]

BC_dissimilarity_components <- ggplot() +
  geom_point(data = dissim_long[`Dissimilarity metric` != "jaccard_dissimilarity_total_mean",], aes(x = year, y = value, color = `Dissimilarity metric`), alpha = 0.2, size = 1) +
  geom_ribbon(data = lmer_bray_total_predictions, aes(x = year, ymin = bray_curtis_lmer_preds_lowerCI, ymax = bray_curtis_lmer_preds_upperCI), fill = "black", alpha = 0.3) +
  geom_line(data = lmer_bray_total_predictions, aes(x = year, y = bray_curtis_lmer_preds), color = "black") +
  #balanced
    geom_ribbon(data = lmer_bray_gradient_predictions, aes(x = year, ymin = bray_curtis_lmer_preds_lowerCI, ymax = bray_curtis_lmer_preds_upperCI), fill = "cadetblue3", alpha = 0.3) +
  geom_line(data = lmer_bray_gradient_predictions, aes(x = year, y = bray_curtis_lmer_preds), color = "cadetblue3") +
  #gradient
    geom_ribbon(data = lmer_bray_balanced_predictions, aes(x = year, ymin = bray_curtis_lmer_preds_lowerCI, ymax = bray_curtis_lmer_preds_upperCI), fill = "darksalmon", alpha = 0.3) +
  geom_line(data = lmer_bray_balanced_predictions, aes(x = year, y = bray_curtis_lmer_preds), color = "darksalmon") +
  #fix. colors of points
  scale_color_manual(values = c("black","darksalmon","cadetblue3"), labels = c("Total","Balanced variation\nin abundance","Biomass gradient")) +
  labs(x = "Year", y = "Value") +
  lims(y = c(min(dissim_long$value), max(dissim_long$value))) +
  theme_classic() +
  theme(legend.position = c(0.17,0.37))

ggsave(BC_dissimilarity_components, path = here::here("figures"), filename=("BC_dissimilarity_components.jpg"), width = 6.5, height = 5, units = "in")

#Jaccard
Jaccard_dissimilarity <- ggplot() +
  geom_point(data = dissim_long[`Dissimilarity metric` == "jaccard_dissimilarity_total_mean"], aes(x = year, y = value), alpha = 0.2, size = 1) +
  geom_ribbon(data = lmer_jaccard_total_predictions, aes(x = year, ymin = jaccard_lmer_preds_lowerCI, ymax = jaccard_lmer_preds_upperCI), fill = "black", alpha = 0.3) +
  geom_line(data = lmer_jaccard_total_predictions, aes(x = year, y = jaccard_lmer_preds), color = "black") +
  labs(x = "Year", y = "Total Jaccard dissimilarity") +
  lims(y = c(min(dissim_long$value), max(dissim_long$value))) +
  theme_classic()

ggsave(Jaccard_dissimilarity, path = here::here("figures"), filename=("Jaccard_dissimilarity.jpg"), width = 6.5, height = 5, units = "in")

#merge
library(cowplot)
all_dissimilarities_over_time <- plot_grid(BC_dissimilarity_components + theme(legend.background = element_rect(fill = "transparent")),
                                           Jaccard_dissimilarity,
                                           ncol = 2, align = "hv", labels = c("a.","b."))

ggsave(all_dissimilarities_over_time, filename = "all_dissimilarities_over_time.jpg", path = here::here("figures"), width = 10, height = 4.5, units = "in")

NA
NA

Bar Plot Coefficient LMER

BC_total_Dissimilarity_Coef_errorbar_reduced <- ggplot() +
    geom_errorbar(data = model_coefs_reduced_length, aes(x = reorder(Survey_Name_Season, year_adj) , y = year_adj, label = Survey_Name_Season, ymin = lwr, ymax = upr), fill = "grey", width = 0) + #add confidence intervals
  geom_point(data = model_coefs_reduced_length, aes(x = reorder(Survey_Name_Season, year_adj) , y = year_adj, label = Survey_Name_Season, size = years_sampled, fill = Directional_Change_sig, color = Directional_Change_sig), stat = 'identity', shape = 21) +
  scale_fill_manual(values = c("white","black","grey"), name = "Dissimilarity trend", guide="none") +
  scale_color_manual(values = c("black","black","grey"), name = "Dissimilarity trend", guide="none") +
  scale_size_binned(range = c(1,8), name = "Survey period length\n(years)") +
  geom_hline(yintercept = 0) +
  scale_y_continuous(breaks = seq(-0.005, 0.0075, by = 0.0025), labels = c("-0.005","","0", "", "0.005",  "")) +
  xlab("Survey unit") +
  ylab("β-diversity trend") + #total BC dissimilarity trend
  coord_flip() +
  theme_classic() +
  theme(axis.text.y = element_text(colour = color_year_adj_order, face = "bold"), axis.title.y = element_blank(), axis.text.x = element_text(size = 15), axis.title.x = element_text(size = 15), legend.position = c(0.25,0.7), legend.direction = "vertical")
Warning: Ignoring unknown parameters: `fill`Warning: Ignoring unknown aesthetics: labelWarning: Ignoring unknown aesthetics: labelWarning: Vectorized input to `element_text()` is not officially supported.
ℹ Results may be unexpected or may change in future versions of ggplot2.
#pull legend for homogenization
directional_change_legend_plot <- BC_total_Dissimilarity_Coef_errorbar_reduced + 
  scale_fill_manual(values = c("white","black","grey"), name = "Dissimilarity trend") +
  scale_color_manual(values = c("black","black","grey"), name = "Dissimilarity trend") +
  scale_size_binned(range = c(1,8), name = "Survey period length", guide = "none") +
  theme(legend.position = "right", legend.background = element_rect(fill= "transparent")) +
  guides(colour = guide_legend(override.aes = list(size=6)))
Scale for fill is already present.
Adding another scale for fill, which will replace the existing scale.Scale for colour is already present.
Adding another scale for colour, which will replace the existing scale.Scale for size is already present.
Adding another scale for size, which will replace the existing scale.
BC_total_Dissimilarity_Coef_errorbar_reduced


#ALT grey scale
BC_total_Dissimilarity_Coef_errorbar_reduced_greyscale <- ggplot() +
    geom_errorbar(data = model_coefs_reduced_length, aes(x = reorder(Survey_Name_Season, year_adj) , y = year_adj, label = Survey_Name_Season, ymin = lwr, ymax = upr), fill = "grey", width = 0) + #add confidence intervals
  geom_point(data = model_coefs_reduced_length, aes(x = reorder(Survey_Name_Season, year_adj) , y = year_adj, label = Survey_Name_Season, size = years_sampled, fill = Directional_Change_sig, color = Directional_Change_sig), stat = 'identity', shape = 21) +
  scale_fill_manual(values = c("white","black","grey"), name = "Dissimilarity trend", guide="none") +
  scale_color_manual(values = c("black","black","grey"), name = "Dissimilarity trend", guide="none") +
  scale_size_binned(range = c(1,8), name = "Survey period length\n(years)") +
  geom_hline(yintercept = 0) +
  scale_y_continuous(breaks = seq(-0.005, 0.0075, by = 0.0025), labels = c("-0.005","","0", "", "0.005",  "")) +
  xlab("Survey unit") +
  ylab("β-diversity trend") + #total BC dissimilarity trend
  coord_flip() +
  theme_classic() +
  theme(axis.text.y = element_text(face = "bold"), axis.title.y = element_blank(), axis.text.x = element_text(size = 15), axis.title.x = element_text(size = 15), legend.position = c(0.25,0.7), legend.direction = "vertical")
Warning: Ignoring unknown parameters: `fill`Warning: Ignoring unknown aesthetics: labelWarning: Ignoring unknown aesthetics: label

Alternatively, we color this plot by trend experienced

#"#73BA4D","#E0962C","#cbbfde"

BC_total_Dissimilarity_Coef_errorbar_reduced_colorbytrend <- ggplot() +
    geom_errorbar(data = model_coefs_reduced_length, aes(x = reorder(Survey_Name_Season, year_adj) , y = year_adj, label = Survey_Name_Season, ymin = lwr, ymax = upr), fill = "grey", width = 0) + #add confidence intervals
  geom_point(data = model_coefs_reduced_length, aes(x = reorder(Survey_Name_Season, year_adj) , y = year_adj, label = Survey_Name_Season, size = years_sampled, color = Directional_Change_sig), stat = 'identity') +
  scale_color_manual(values = c("#73BA4D","#E0962C","#cbbfde"), name = "Dissimilarity trend", guide="none") +
  scale_size_binned(range = c(1,8), name = "Survey period length\n(years)") +
  geom_hline(yintercept = 0) +
  scale_y_continuous(breaks = seq(-0.005, 0.0075, by = 0.0025), labels = c("-0.005","","0", "", "0.005",  "")) +
  xlab("Survey unit") +
  ylab("β-diversity trend") + #total BC dissimilarity trend
  coord_flip() +
  theme_classic() +
  theme(axis.text.y = element_text(face = "bold"), axis.title.y = element_blank(), axis.text.x = element_text(size = 15), axis.title.x = element_text(size = 15), legend.position = c(0.3,0.8), legend.direction = "vertical", legend.text = element_text(size = 15), legend.title = element_text(size = 16))
Warning: Ignoring unknown parameters: `fill`Warning: Ignoring unknown aesthetics: labelWarning: Ignoring unknown aesthetics: label
directional_change_legend_plot_colorbytrend <- BC_total_Dissimilarity_Coef_errorbar_reduced_colorbytrend + 
  theme(legend.position = "right", legend.background = element_rect(fill= "transparent"), 
         legend.text = element_text(size = 15), legend.title = element_text(size = 16)) +
  guides(colour = guide_legend(override.aes = list(size=6)), size = "none")

Wavy Line Plot for GAMs

Generate predicted values


#add colors and names to full dissimilarity data table
distances_dissimilarities_allyears.r <- distances_dissimilarities_allyears.r[color_link, on = "survey_unit"]

#generate new data to smooth lines (need year and season survey combinations)
year_survey_unit_expand.dt <- data.table(survey_unit = as.character(NULL), year = as.numeric(NULL), year_adj = as.numeric(NULL ))

for (i in 1:length(survey_unit.list)) {
  #generate year vectors
  survey_unit_years <- unique(distances_dissimilarities_allyears.r[survey_unit == survey_unit.list[i],.(survey_unit, year, year_adj)])
  
  years <- seq(min(survey_unit_years[,year]), max(survey_unit_years[,year]), by = 0.1)
  
  year_adjust <- seq(min(survey_unit_years[,year_adj]), max(survey_unit_years[,year_adj]), by = 0.1)
  
  year_survey_unit_expand.dt_addition <- data.table(survey_unit = survey_unit.list[i], year = years, year_adj = year_adjust)
  
  year_survey_unit_expand.dt <- rbind(year_survey_unit_expand.dt, year_survey_unit_expand.dt_addition)
}

#add colors and names to full year and survey unit combination table
year_survey_unit_expand.dt <- year_survey_unit_expand.dt[color_link, on = "survey_unit"]

Get model as predictions

#for plotting, get model as predictions
bray_curtis_total_gam_predictions <- predict(bray_curtis_total_gam, se.fit = TRUE, newdata = year_survey_unit_expand.dt)

#merge into table
year_survey_unit_expand.dt[,bray_glm_mod_fit := bray_curtis_total_gam_predictions$fit][,bray_glm_mod_fit_SE := bray_curtis_total_gam_predictions$se.fit]

Produce Plot of GAM and mean LMER line

points_wavylines_bray_total_year_reduced_gam_nolmer <- ggplot() +
 # geom_ribbon(data = lmer_bray_total_predictions, aes(x = year, ymin = bray_curtis_lmer_preds_lowerCI, ymax = bray_curtis_lmer_preds_upperCI), fill = "grey", alpha = 0.2) +
  geom_point(data = na.omit(distances_dissimilarities_allyears.r,cols = "year_adj"),
             aes(x = year,
                 y = bray_curtis_dissimilarity_total_mean,
                 color = Survey_Name_Season), alpha = 0.5, size = 1) +
    geom_line(data = na.omit(year_survey_unit_expand.dt,cols = "year_adj"),
             aes(x = year,
                 y = bray_glm_mod_fit,
                 color = Survey_Name_Season)) +
  geom_ribbon(data = na.omit(year_survey_unit_expand.dt,cols = "year_adj"), aes(x = year, ymin=bray_glm_mod_fit-bray_glm_mod_fit_SE, ymax=bray_glm_mod_fit+bray_glm_mod_fit_SE, fill =  Survey_Name_Season), alpha=0.1) + #add standard error
 # geom_line(data = lmer_bray_total_predictions, aes(x = year, y = bray_curtis_lmer_preds), color = "black") +
    scale_color_manual(values =  color_alpha_order, name = "Survey Unit") +
  scale_fill_manual(values =  color_alpha_order, guide = "none") +
  theme_classic() +
  lims(x = c(min(distances_dissimilarities_allyears.r[,year]),max(distances_dissimilarities_allyears.r[,year])),
       y = c(0.15,0.9)) +
  xlab("Year") +
ylab("total BC dissimilarity") +
  theme(legend.position = "null")

points_wavylines_bray_total_year_reduced_gam_nolmer

ggsave(points_wavylines_bray_total_year_reduced_gam_nolmer, path = here::here("figures"), filename ="points_wavylines_bray_total_year_reduced_gam_nolmer.jpg", height = 5, width = 5, unit = "in")


#with lmer

points_wavylines_bray_total_year_reduced_gam <- ggplot() +
  geom_ribbon(data = lmer_bray_total_predictions, aes(x = year, ymin = bray_curtis_lmer_preds_lowerCI, ymax = bray_curtis_lmer_preds_upperCI), fill = "grey", alpha = 0.3) +
  geom_point(data = na.omit(distances_dissimilarities_allyears.r, cols = "year_adj"),
             aes(x = year,
                 y = bray_curtis_dissimilarity_total_mean,
                 fill = Survey_Name_Season), alpha = 0.4, size = 1, shape = 21, color = "white") +
    geom_line(data = na.omit(year_survey_unit_expand.dt, cols = "year_adj"),
             aes(x = year,
                 y = bray_glm_mod_fit,
                 color = Survey_Name_Season), alpha = 0.6) +
  geom_ribbon(data = na.omit(year_survey_unit_expand.dt, cols = "year_adj"), aes(x = year, ymin=bray_glm_mod_fit-bray_glm_mod_fit_SE, ymax=bray_glm_mod_fit+bray_glm_mod_fit_SE, fill =  Survey_Name_Season), alpha=0.1) + #add standard error
  geom_line(data = lmer_bray_total_predictions, aes(x = year, y = bray_curtis_lmer_preds), color = "black") +
    scale_color_manual(values =  color_alpha_order, name = "Survey Unit") +
  scale_fill_manual(values =  color_alpha_order, guide = "none") +
  theme_classic() +
  lims(x = c(min(distances_dissimilarities_allyears.r[,year]),max(distances_dissimilarities_allyears.r[,year]))) +
  xlab("Year") +
ylab("β-diversity") +
  theme(legend.position = "null", axis.text = element_text(size = 15), axis.title = element_text(size = 15))

points_wavylines_bray_total_year_reduced_gam

ggsave(points_wavylines_bray_total_year_reduced_gam, path = here::here("figures"), filename ="points_wavylines_bray_total_year_reduced_gam.jpg", height = 6, width = 6, unit = "in")


#ALT
#plot all, but same color scheme (grey)
points_wavylines_bray_total_year_reduced_gam_greyscale <- ggplot() +
  geom_ribbon(data = lmer_bray_total_predictions, aes(x = year, ymin = bray_curtis_lmer_preds_lowerCI, ymax = bray_curtis_lmer_preds_upperCI), fill = "grey", alpha = 0.3) +
  geom_point(data = distances_dissimilarities_allyears.r,
             aes(x = year,
                 y = bray_curtis_dissimilarity_total_mean,
                 fill = Survey_Name_Season), alpha = 0.4, size = 1, shape = 21, color = "white") +
    geom_line(data = year_survey_unit_expand.dt,
             aes(x = year,
                 y = bray_glm_mod_fit,
                 color = Survey_Name_Season), alpha = 0.6) +
  geom_ribbon(data = year_survey_unit_expand.dt, aes(x = year, ymin=bray_glm_mod_fit-bray_glm_mod_fit_SE, ymax=bray_glm_mod_fit+bray_glm_mod_fit_SE, fill =  Survey_Name_Season), alpha=0.1) + #add standard error
  geom_line(data = lmer_bray_total_predictions, aes(x = year, y = bray_curtis_lmer_preds), color = "black") +
    scale_color_manual(values =  rep("black", times = length(unique(distances_dissimilarities_allyears.r$Survey_Name_Season))), name = "Survey Unit") +
  scale_fill_manual(values =  rep("black", times = length(unique(distances_dissimilarities_allyears.r$Survey_Name_Season))), guide = "none") +
  theme_classic() +
  lims(x = c(min(distances_dissimilarities_allyears.r[,year]),max(distances_dissimilarities_allyears.r[,year]))
       ) +
  xlab("Year") +
ylab("β-diversity") +
  theme(legend.position = "null", axis.text = element_text(size = 15), axis.title = element_text(size = 15))

points_wavylines_bray_total_year_reduced_gam_greyscale

ggsave(points_wavylines_bray_total_year_reduced_gam_greyscale, path = here::here("figures"), filename ="points_wavylines_bray_total_year_reduced_gam_greyscale.jpg", height = 6, width = 6, unit = "in")

Alternative, color by trend


points_wavylines_bray_total_year_reduced_gam_colorbytrend <- ggplot() +
  geom_ribbon(data = lmer_bray_total_predictions, aes(x = year, ymin = bray_curtis_lmer_preds_lowerCI, ymax = bray_curtis_lmer_preds_upperCI), fill = "grey", alpha = 0.3) +
  geom_point(data = na.omit(distances_dissimilarities_allyears.r, cols = "year_adj"),
             aes(x = year,
                 y = bray_curtis_dissimilarity_total_mean,
                 fill = Survey_Name_Season), alpha = 0.4, size = 1, shape = 21, color = "white") +
    geom_line(data = na.omit(year_survey_unit_expand.dt, cols = "year_adj"),
             aes(x = year,
                 y = bray_glm_mod_fit,
                 color = Survey_Name_Season), alpha = 0.6) +
  geom_ribbon(data = na.omit(year_survey_unit_expand.dt, cols = "year_adj"), aes(x = year, ymin=bray_glm_mod_fit-bray_glm_mod_fit_SE, ymax=bray_glm_mod_fit+bray_glm_mod_fit_SE, fill =  Survey_Name_Season), alpha=0.1) + #add standard error
  geom_line(data = lmer_bray_total_predictions, aes(x = year, y = bray_curtis_lmer_preds), color = "black") +
    scale_color_manual(values =  color_alpha_order_bytrend, name = "Survey Unit") +
  scale_fill_manual(values =  color_alpha_order_bytrend, guide = "none") +
  theme_classic() +
  lims(x = c(min(distances_dissimilarities_allyears.r[,year]),max(distances_dissimilarities_allyears.r[,year])),y = c(0.49,0.95)) +
  xlab("Year") +
ylab("β-diversity") +
  theme(legend.position = "null", axis.text = element_text(size = 15), axis.title = element_text(size = 15))

points_wavylines_bray_total_year_reduced_gam_colorbytrend

ggsave(points_wavylines_bray_total_year_reduced_gam_colorbytrend, path = here::here("figures"), filename ="points_wavylines_bray_total_year_reduced_gam_colorbytrend.jpg", height = 6, width = 11, unit = "in")



points_wavylines_bray_total_year_reduced_gam_colorbytrend_squigglesonly <- ggplot() +
  geom_point(data = na.omit(distances_dissimilarities_allyears.r, cols = "year_adj"),
             aes(x = year,
                 y = bray_curtis_dissimilarity_total_mean,
                 fill = Survey_Name_Season), alpha = 0.4, size = 1, shape = 21, color = "white") +
    geom_line(data = na.omit(year_survey_unit_expand.dt, cols = "year_adj"),
             aes(x = year,
                 y = bray_glm_mod_fit,
                 color = Survey_Name_Season), alpha = 0.6) +
  geom_ribbon(data = na.omit(year_survey_unit_expand.dt, cols = "year_adj"), aes(x = year, ymin=bray_glm_mod_fit-bray_glm_mod_fit_SE, ymax=bray_glm_mod_fit+bray_glm_mod_fit_SE, fill =  Survey_Name_Season), alpha=0.1) + #add standard error
    scale_color_manual(values =  color_alpha_order_bytrend, name = "Survey Unit") +
  scale_fill_manual(values =  color_alpha_order_bytrend, guide = "none") +
  theme_classic() +
  lims(x = c(min(distances_dissimilarities_allyears.r[,year]),max(distances_dissimilarities_allyears.r[,year])), y = c(0.49,0.95)) +
  xlab("Year") +
ylab("β-diversity") +
  theme(legend.position = "null", axis.text = element_text(size = 15), axis.title = element_text(size = 15))

points_wavylines_bray_total_year_reduced_gam_colorbytrend_squigglesonly

ggsave(points_wavylines_bray_total_year_reduced_gam_colorbytrend_squigglesonly, path = here::here("figures"), filename ="points_wavylines_bray_total_year_reduced_gam_colorbytrend_squigglesonly.jpg", height = 6, width = 11, unit = "in")

#plot each independently for supplement
#all survey names = 
all_survey_names <- sort(unique(color_link$Survey_Name_Season))
#list of plots
points_wavylines_bray_total_year_reduced_gam_individual <- list()
for (i in 1:length(all_survey_names)) {
points_wavylines_bray_total_year_reduced_gam_individual[[i]] <- ggplot() +
  geom_point(data = distances_dissimilarities_allyears.r[Survey_Name_Season == all_survey_names[i]],
             aes(x = year,
                 y = bray_curtis_dissimilarity_total_mean), alpha = 0.4, color = "black") +
    geom_line(data = year_survey_unit_expand.dt[Survey_Name_Season == all_survey_names[i]],
             aes(x = year,
                 y = bray_glm_mod_fit), alpha = 0.6) +
  geom_ribbon(data = year_survey_unit_expand.dt[Survey_Name_Season == all_survey_names[i]], aes(x = year, ymin=bray_glm_mod_fit-bray_glm_mod_fit_SE, ymax=bray_glm_mod_fit+bray_glm_mod_fit_SE), alpha=0.1) + #add standard error
  theme_classic() +
#  lims(x = c(min(distances_dissimilarities_allyears.r[Survey_Name_Season == all_survey_names[i],year]),max(distances_dissimilarities_allyears.r[Survey_Name_Season == all_survey_names[i],year])),
#       y = c(0.15,0.9)) +
  xlab("Year") +
ylab("beta-diversity") +
  facet_wrap(~Survey_Name_Season, ncol = 5) +
  theme(legend.position = "null", axis.text = element_text(size = 15), axis.title = element_text(size = 15))

print(points_wavylines_bray_total_year_reduced_gam_individual[[i]])

}

saveRDS(points_wavylines_bray_total_year_reduced_gam_individual, here::here("figures","points_wavylines_bray_total_year_reduced_gam_individual.Rds"))

#print to pdf
library(gridExtra)

ggsave(
   filename = here::here("figures","points_wavylines_bray_total_year_reduced_gam_individual.pdf"), 
   plot = marrangeGrob(points_wavylines_bray_total_year_reduced_gam_individual, nrow=1, ncol=1), 
   width = 8.5, height = 11
)

#Alternatively, split into 2 and use facet
#first 24
points_wavylines_bray_total_year_reduced_gam_individual_facet_1_24 <- ggplot() +
  geom_point(data = distances_dissimilarities_allyears.r[Survey_Name_Season  %in% all_survey_names[c(1:24)]],
             aes(x = year,
                 y = bray_curtis_dissimilarity_total_mean), alpha = 0.7) +
    geom_line(data = year_survey_unit_expand.dt[Survey_Name_Season   %in% all_survey_names[c(1:24)]],
             aes(x = year,
                 y = bray_glm_mod_fit)) +
  geom_ribbon(data = year_survey_unit_expand.dt[Survey_Name_Season  %in% all_survey_names[c(1:24)]],
aes(x = year, ymin=bray_glm_mod_fit-bray_glm_mod_fit_SE, ymax=bray_glm_mod_fit+bray_glm_mod_fit_SE), alpha=0.5) + #add standard error
  theme_classic() +
  xlab("Year") +
ylab("β-diversity") +
  facet_wrap(~Survey_Name_Season, ncol = 4, scales = "free") +
  theme(axis.text = element_text(size = 8), axis.title = element_text(size = 12))

ggsave(points_wavylines_bray_total_year_reduced_gam_individual_facet_1_24, path =  here::here("figures"), filename = "points_wavylines_bray_total_year_reduced_gam_individual_facet_1_24.png", height = 11, width = 9)

points_wavylines_bray_total_year_reduced_gam_individual_facet_25_42 <- ggplot() +
  geom_point(data = distances_dissimilarities_allyears.r[Survey_Name_Season  %in% all_survey_names[c(25:42)]],
             aes(x = year,
                 y = bray_curtis_dissimilarity_total_mean), alpha = 0.7) +
    geom_line(data = year_survey_unit_expand.dt[Survey_Name_Season   %in% all_survey_names[c(25:42)]],
             aes(x = year,
                 y = bray_glm_mod_fit)) +
  geom_ribbon(data = year_survey_unit_expand.dt[Survey_Name_Season  %in% all_survey_names[c(25:42)]],
aes(x = year, ymin=bray_glm_mod_fit-bray_glm_mod_fit_SE, ymax=bray_glm_mod_fit+bray_glm_mod_fit_SE), alpha=0.5) + #add standard error
  theme_classic() +
  xlab("Year") +
ylab("β-diversity") +
  facet_wrap(~Survey_Name_Season, ncol = 4, scales = "free") +
  theme(axis.text = element_text(size = 8), axis.title = element_text(size = 12))

ggsave(points_wavylines_bray_total_year_reduced_gam_individual_facet_25_42, path =  here::here("figures"), filename = "points_wavylines_bray_total_year_reduced_gam_individual_facet_25_42.png", height = 11, width = 9)

NA
NA

Region by region build up

#EBS ONLY (segue for chapter 3)
points_wavylines_bray_total_year_reduced_gam_colorbytrend_EBSonly <- ggplot() +
  geom_point(data = na.omit(distances_dissimilarities_allyears.r[survey_unit == "EBS"], cols = "year_adj"),
             aes(x = year,
                 y = bray_curtis_dissimilarity_total_mean,
                 fill = Survey_Name_Season), size = 2, shape = 21, color = "white") +
    geom_line(data = na.omit(year_survey_unit_expand.dt[survey_unit == "EBS"], cols = "year_adj"),
             aes(x = year,
                 y = bray_glm_mod_fit,
                 color = Survey_Name_Season)) +
  geom_ribbon(data = na.omit(year_survey_unit_expand.dt[survey_unit == "EBS"], cols = "year_adj"), aes(x = year, ymin=bray_glm_mod_fit-bray_glm_mod_fit_SE, ymax=bray_glm_mod_fit+bray_glm_mod_fit_SE, fill =  Survey_Name_Season), alpha=0.1) + #add standard error
    scale_color_manual(values =  "#cbbfde", name = "Survey Unit") +
  scale_fill_manual(values =  "#cbbfde", guide = "none") +
  theme_classic() +
  lims(x = c(min(distances_dissimilarities_allyears.r[,year]),max(distances_dissimilarities_allyears.r[,year])), y = c(0.49,0.95)) +
  xlab("Year") +
ylab("β-diversity") +
  theme(legend.position = "null", axis.text = element_text(size = 15), axis.title = element_text(size = 15))

points_wavylines_bray_total_year_reduced_gam_colorbytrend_EBSonly

ggsave(points_wavylines_bray_total_year_reduced_gam_colorbytrend_EBSonly, path = here::here("figures"), filename ="points_wavylines_bray_total_year_reduced_gam_colorbytrend_EBSonly.jpg", height = 6, width = 11, unit = "in")



#just one region that's homogenizing (SEUS-summer)
points_wavylines_bray_total_year_reduced_gam_colorbytrend_seus_summer <- ggplot() +
  geom_point(data = na.omit(distances_dissimilarities_allyears.r[survey_unit == "SEUS-summer"], cols = "year_adj"),
             aes(x = year,
                 y = bray_curtis_dissimilarity_total_mean,
                 fill = Survey_Name_Season), size = 2, shape = 21, color = "white") +
    geom_line(data = na.omit(year_survey_unit_expand.dt[survey_unit == "SEUS-summer"], cols = "year_adj"),
             aes(x = year,
                 y = bray_glm_mod_fit,
                 color = Survey_Name_Season)) +
  geom_ribbon(data = na.omit(year_survey_unit_expand.dt[survey_unit == "SEUS-summer"], cols = "year_adj"), aes(x = year, ymin=bray_glm_mod_fit-bray_glm_mod_fit_SE, ymax=bray_glm_mod_fit+bray_glm_mod_fit_SE, fill =  Survey_Name_Season), alpha=0.1) + #add standard error
    scale_color_manual(values =  "#e7ac5b", name = "Survey Unit") +
  scale_fill_manual(values =  "#e7ac5b", guide = "none") +
  theme_classic() +
  lims(x = c(min(distances_dissimilarities_allyears.r[,year]),max(distances_dissimilarities_allyears.r[,year])), y = c(0.49,0.95)) +
  xlab("Year") +
ylab("β-diversity") +
  theme(legend.position = "null", axis.text = element_text(size = 15), axis.title = element_text(size = 15))

points_wavylines_bray_total_year_reduced_gam_colorbytrend_seus_summer

ggsave(points_wavylines_bray_total_year_reduced_gam_colorbytrend_seus_summer, path = here::here("figures"), filename ="points_wavylines_bray_total_year_reduced_gam_colorbytrend_seus_summer.jpg", height = 6, width = 11, unit = "in")




#just one region that's differentiating and one that's homogenizing
points_wavylines_bray_total_year_reduced_gam_colorbytrend_NZ_WCSI_SEUS_summeronly <- ggplot() +
  geom_point(data = na.omit(distances_dissimilarities_allyears.r[survey_unit %in% c("NZ-WCSI","SEUS-summer")], cols = "year_adj"),
             aes(x = year,
                 y = bray_curtis_dissimilarity_total_mean,
                 fill = Survey_Name_Season), size = 2, shape = 21, color = "white") +
    geom_line(data = na.omit(year_survey_unit_expand.dt[survey_unit %in% c("NZ-WCSI","SEUS-summer")], cols = "year_adj"),
             aes(x = year,
                 y = bray_glm_mod_fit,
                 color = Survey_Name_Season)) +
  geom_ribbon(data = na.omit(year_survey_unit_expand.dt[survey_unit %in% c("NZ-WCSI","SEUS-summer")], cols = "year_adj"), aes(x = year, ymin=bray_glm_mod_fit-bray_glm_mod_fit_SE, ymax=bray_glm_mod_fit+bray_glm_mod_fit_SE, fill =  Survey_Name_Season), alpha=0.1) + #add standard error
    scale_color_manual(values =  c("#e7ac5b","#91c874"), name = "Survey Unit") +
  scale_fill_manual(values =  c("#e7ac5b","#91c874"), guide = "none") +
  theme_classic() +
  lims(x = c(min(distances_dissimilarities_allyears.r[,year]),max(distances_dissimilarities_allyears.r[,year])), y = c(0.49,0.95)) +
  xlab("Year") +
ylab("β-diversity") +
  theme(legend.position = "null", axis.text = element_text(size = 15), axis.title = element_text(size = 15))

points_wavylines_bray_total_year_reduced_gam_colorbytrend_NZ_WCSI_SEUS_summeronly

ggsave(points_wavylines_bray_total_year_reduced_gam_colorbytrend_NZ_WCSI_SEUS_summeronly, path = here::here("figures"), filename ="points_wavylines_bray_total_year_reduced_gam_colorbytrend_NZ_WCSI_SEUS_summeronly.jpg", height = 6, width = 11, unit = "in")


#just three regions, another one that's stable, (NEUS-Fall)
points_wavylines_bray_total_year_reduced_gam_colorbytrend_NZ_WCSI_SEUS_summer_NEUSfallonly <- ggplot() +
  geom_point(data = na.omit(distances_dissimilarities_allyears.r[survey_unit %in% c("NZ-WCSI","SEUS-summer", "NEUS-Fall")], cols = "year_adj"),
             aes(x = year,
                 y = bray_curtis_dissimilarity_total_mean,
                 fill = Survey_Name_Season), size = 2, shape = 21, color = "white") +
    geom_line(data = na.omit(year_survey_unit_expand.dt[survey_unit %in% c("NZ-WCSI","SEUS-summer", "NEUS-Fall")], cols = "year_adj"),
             aes(x = year,
                 y = bray_glm_mod_fit,
                 color = Survey_Name_Season)) +
  geom_ribbon(data = na.omit(year_survey_unit_expand.dt[survey_unit %in% c("NZ-WCSI","SEUS-summer", "NEUS-Fall")], cols = "year_adj"), aes(x = year, ymin=bray_glm_mod_fit-bray_glm_mod_fit_SE, ymax=bray_glm_mod_fit+bray_glm_mod_fit_SE, fill =  Survey_Name_Season), alpha=0.1) + #add standard error
    scale_color_manual(values =  c("#cbbfde","#e7ac5b", "#91c874"), name = "Survey Unit") +
  scale_fill_manual(values =  c("#cbbfde","#e7ac5b", "#91c874"), guide = "none") +
  theme_classic() +
  lims(x = c(min(distances_dissimilarities_allyears.r[,year]),max(distances_dissimilarities_allyears.r[,year])), y = c(0.49,0.95)) +
  xlab("Year") +
ylab("β-diversity") +
  theme(legend.position = "null", axis.text = element_text(size = 15), axis.title = element_text(size = 15))

points_wavylines_bray_total_year_reduced_gam_colorbytrend_NZ_WCSI_SEUS_summer_NEUSfallonly

ggsave(points_wavylines_bray_total_year_reduced_gam_colorbytrend_NZ_WCSI_SEUS_summer_NEUSfallonly, path = here::here("figures"), filename ="points_wavylines_bray_total_year_reduced_gam_colorbytrend_NZ_WCSI_SEUS_summer_NEUSfallonly.jpg", height = 6, width = 11, unit = "in")

NA
NA
NA

Merge BC versus Year plot with GAMS and Region vs. coefficient plot for LMERs


BC_total_GAM_LMER_merge_legend <- ggdraw(xlim = c(0, 40.5), ylim = c(0, 21)) +
    draw_plot(points_wavylines_bray_total_year_reduced_gam,
                                         x = 1, y = 1, width = 20, height = 20) +
    draw_plot(BC_total_Dissimilarity_Coef_errorbar_reduced +
        theme(legend.key.size = unit(0.5, 'cm'), #change legend key size
        legend.title = element_text(size=16), #change legend title font size
        legend.text = element_text(size=14)), #change legend text font size),
                                         x = 20, y = 1, width =19, height = 20) +
    draw_plot(get_legend(directional_change_legend_plot + 
      theme(legend.key.size = unit(0.5, 'cm'), #change legend key size
        legend.title = element_text(size=15), #change legend title font size
        legend.text = element_text(size=13))), #change legend text font size)
                                       x = 27, y = 10, width = 3, height = 2) +
  geom_text(aes(x = 2, y = 20.7), label = ("a."), size =8, fontface = "bold") +
  geom_text(aes(x =20, y = 20.7), label = ("b."), size =8, fontface = "bold")


ggsave(BC_total_GAM_LMER_merge_legend, path = here::here("figures"), filename = "BC_total_GAM_LMER_merge_legend.png", height = 8, width = 14, units = "in")

#ALT GREY SCALE
BC_total_GAM_LMER_merge_legend_greyscale <- ggdraw(xlim = c(0, 40.5), ylim = c(0, 21)) +
    draw_plot(points_wavylines_bray_total_year_reduced_gam_greyscale,
                                         x = 1, y = 1, width = 20, height = 20) +
    draw_plot(BC_total_Dissimilarity_Coef_errorbar_reduced_greyscale +
        theme(legend.key.size = unit(0.5, 'cm'), #change legend key size
        legend.title = element_text(size=16), #change legend title font size
        legend.text = element_text(size=14)), #change legend text font size),
                                         x = 20, y = 1, width = 19, height = 20) +
    draw_plot(get_legend(directional_change_legend_plot + 
      theme(legend.key.size = unit(0.5, 'cm'), #change legend key size
        legend.title = element_text(size=15), #change legend title font size
        legend.text = element_text(size=13))), #change legend text font size)
                                x = 27, y = 10, width = 3, height = 2) +
  geom_text(aes(x = 2, y = 20.7), label = ("a."), size =8, fontface = "bold") +
  geom_text(aes(x =20, y = 20.7), label = ("b."), size =8, fontface = "bold")

ggsave(BC_total_GAM_LMER_merge_legend_greyscale, path = here::here("figures"), filename = "BC_total_GAM_LMER_merge_legend_greyscale.png", height = 8, width = 14, units = "in")

#ALT COLOR BY TREND
BC_total_GAM_LMER_merge_legend_colorbytrend <- ggdraw(xlim = c(0, 40.5), ylim = c(0, 21)) +
    draw_plot(points_wavylines_bray_total_year_reduced_gam_colorbytrend,
                                         x = 1, y = 1, width = 20, height = 20) +
    draw_plot(BC_total_Dissimilarity_Coef_errorbar_reduced_colorbytrend +
        theme(legend.key.size = unit(0.5, 'cm'), #change legend key size
       # legend.title = element_text(size=16), #change legend title font size
       # legend.text = element_text(size=14)
       ), #change legend text font size),
                                         x = 20, y = 1, width = 19, height = 20) +
    draw_plot(get_legend(directional_change_legend_plot_colorbytrend + 
      theme(legend.key.size = unit(0.5, 'cm'), #change legend key size
        legend.title = element_text(size=16), #change legend title font size
        legend.text = element_text(size=15))), #change legend text font size)
                                x = 27, y = 8, width = 3, height = 2) +
  geom_text(aes(x = 2, y = 20.7), label = ("a."), size =8, fontface = "bold") +
  geom_text(aes(x =20, y = 20.7), label = ("b."), size =8, fontface = "bold")

ggsave(BC_total_GAM_LMER_merge_legend_colorbytrend, path = here::here("figures"), filename = "BC_total_GAM_LMER_merge_legend_colorbytrend.png", height = 8, width = 14, units = "in")

Is each survey better described by lmer or by GAM?

LS0tCnRpdGxlOiAiWWVhciBUT1RBTCBEaXNzaW1pbGFyaXR5IE1vZGVscyIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKVGhpcyBjb2RlIGlzIFNjcmlwdCA0YSAgZm9yIEtpdGNoZWwgZXQgYWwuIEJpb3RpYyBob21vZ2VuaXphdGlvbiwgdGhlIGV4Y2VwdGlvbiBhbmQgbm90IHRoZSBydWxlIGZvciBtYXJpbmUgZmlzaCBjb21tdW5pdGllcyBtYW51c2NyaXB0LgoKQ09QWSBGUk9NIFRISVMgU0NSSVBUCgooV2UgZXhwZWN0IHRoaXMgdG8gYmUgdGhlIG1haW4gYW5hbHlzZXMgaW4gdGhlIHRleHQpCgotIFRoaXMgcHJvamVjdCBpcyBhIGNvbGxhYm9yYXRpdmUgZWZmb3J0IHRvIGRlc2NyaWJlIGNoYW5nZXMgaW4gdGF4b25vbWljIGNvbXBvc2l0aW9uICBvZiBmaXNoIGNvbW11bml0aWVzIGFyb3VuZCB0aGUgd29ybGQtLWFzIHNhbXBsZWQgYnkgYm90dG9tIHRyYXdsIHN1cnZleXMuCgotIENvZGUgYnkgWm/DqyBKLiBLaXRjaGVsCgpTRVNTSU9OIElORk8gVE8gRE8KCmBgYHtyIHNldHVwfQpsaWJyYXJ5KGRhdGEudGFibGUpCmxpYnJhcnkodmVnYW4pCmxpYnJhcnkoc2YpCmxpYnJhcnkoY29uY2F2ZW1hbikgI3BvbHlnb24gYXJvdW5kIHBvaW50cwpsaWJyYXJ5KGJldGFwYXJ0KSAjYWxsb3dzIHVzIHRvIHBhcnRpdGlvbiBiZXRhIGRpdmVyc2l0eQpsaWJyYXJ5KGdlb3NwaGVyZSkKbGlicmFyeShnZ3B1YnIpICNzdGF0X3JlZ2xpbmVfZXF1YXRpb24KbGlicmFyeShubG1lKQpsaWJyYXJ5KG1nY3YpICN0byBtYWtlIGdhbQpsaWJyYXJ5KGNvd3Bsb3QpCmxpYnJhcnkobG1lNCkKCiNQdWxsIERpc3NpbWlsYXJpdHkgTWVhbnMKZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yIDwtIHJlYWRSRFMoaGVyZTo6aGVyZSgib3V0cHV0IiwgImRpc3NpbWlsYXJpdGllcyIsICJkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnIucmRzIikpCgojbWFrZSBzdXJ2ZXkgYW5kIHN1cnZleSB1bml0IGZhY3RvcnMKZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyxzdXJ2ZXk6PWZhY3RvcihzdXJ2ZXkpXVssc3VydmV5X3VuaXQ6PWZhY3RvcihzdXJ2ZXlfdW5pdCldCgojYWRqdXN0IHllYXJzCmRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclsseWVhcl9hZGogOj0geWVhci1taW4oeWVhcikrMV0KCiNhZGQgbmV3IHZhcmlhYmxlIGZvciB5ZWFyIGluIHNlcXVlbmNlIHBlciByZWdpb24KZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyxmaXJzdF95ZWFyIDo9IG1pbih5ZWFyKSwuKHN1cnZleV91bml0KV0KZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyxsYXN0X3llYXIgOj0gbWF4KHllYXIpLC4oc3VydmV5X3VuaXQpXQoKI2Rpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclsseWVhcl9pbl9zZXEgOj0geWVhci1maXJzdF95ZWFyKzFdCgpkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLHllYXJzX3NhbXBsZWQgOj0gbGFzdF95ZWFyLWZpcnN0X3llYXIrMV0KCgpgYGAKCiMjI1BhbGV0dGUgZm9yIFBsb3R0aW5nClBhbGV0dGUgZm9yIHBsb3R0aW5nIGFsbCA0MiBzdXJ2ZXkgdW5pdHMKYGBge3IgbGluayBjb2xvcnMgdG8gc3VydmV5IHVuaXRzfQpzdXJ2ZXlfdW5pdC5saXN0IDwtIGxldmVscyhmYWN0b3IoZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyxzdXJ2ZXlfdW5pdF0pKQoKcGFsZXR0ZV80MiA8LSBjKAogICIjNUE1MTU2IiwgI0FJCiAgIiNERjAwREIiLCAjQklUUy0xCiAgIiNEQjhFREEiLCAjQklUUy00CiAgIiNGNjIyMkUiLCAjQ0hMCiAgIiNGOEExOUYiLCAjREZPLU5GCiAgIiMxNkZGMzIiLCAjREZPLVFDUwogICIjMzI1QTlCIiwgI0VCUwogICIjMzI4M0ZFIiwgI0VWSE9FCiAgIiNGRUFGMTYiLCAjRlItQ0dGUwogICIjZmNjYjZkIiwgI0dNRVgtRmFsbAogICIjMUM4MzU2IiwgI0dNRVgtU3VtbWVyCiAgIiNDNDQ1MUMiLCAjR09BCiAgIiM4NTY2MEQiLCAjR1JMLURFCiAgIiNCMDAwOUYiLCAjR1NMLU4KICAiI0JGNzlCOCIsICNHU0wtUwogICIjMUNCRTRGIiwgI0lDRS1HRlMKICAiIzc4MkFCNiIsICNJRS1JR0ZTCiAgIiM5MEFEMUMiLCAjTUVESVRTCiAgIiM2QjAwM0EiLCAjTkFNCiAgIiNBNzVCMDAiLCAjTkVVUy1GYWxsCiAgIiNFM0IwNzIiLCAjTkVVUy1TcHJpbmcKICAiIzAyRThCNiIsICNOSUdGUy0xCiAgIiM5N0U3RDUiLCAjTklHRlMtNAogICIjQjAwMDY4IiwgI05vci1CVFMtMwogICIjMDBCOUUzIiwgI05TLUlCVFMtMQogICIjOTVFMkY0IiwgI05TLUlCVFMtMwogICIjQjNDRTczIiwgI05aLUNIQVQKICAiIzY4OTUwMCIsICNOWi1FQ1NJCiAgIiMzNjRkMDIiLCNOWi1TVUJBCiAgIiNBQUY0MDAiLCAjTlotV0NTSQogICIjQUEwREZFIiwgI1BULUlCVFMKICAiIzdmOWViOCIsICNST0NLQUxMCiAgIiNGQTAwODciLCAjUy1HRU9SRwogICIjREVBMEZEIiwgI1NDUy1TdW1tZXIKICAiI0ZDRUY4OCIsICNTRVVTLWZhbGwKICAiI0E1OTQwNSIsICNTRVVTLXNwcmluZwogICIjRkNFMTAwIiwgI1NFVVMtc3VtbWVyCiAgIiM1NDQ1NjMiLCAjU1dDLUlCVFMtMQogICIjYTM3ZmM3IiwgI1NXQy1JQlRTLTQKICAiI0MwNzVBNiIsICNXQ0FOTgogICIjQkRDREZGIiwgI1pBRi1BVEwKICAiIzAwM0VGRiIgICNaQUYtSU5ECikKCmNvbG9yX2xpbmsgPC0gZGF0YS50YWJsZShzdXJ2ZXlfdW5pdCA9IHN1cnZleV91bml0Lmxpc3QsaGV4ID0gcGFsZXR0ZV80MikKYGBgCgpBZGQgbmFtZXMgZm9yIHBsb3R0aW5nCmBgYHtyIGFkZCBuYW1lcyBmb3IgcGxvdHRpbmd9CgpuYW1lX2hlbHBlciA8LSBkYXRhLnRhYmxlKFN1cnZleV9OYW1lX1NlYXNvbiA9IGMoIkFsZXV0aWFuIElzbGFuZHMiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiQmFsdGljIFNlYSBRMSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJCYWx0aWMgU2VhIFE0IiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkNoaWxlIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5ld2ZvdW5kbGFuZCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJRdWVlbiBDaGFybG90dGUgU291bmQiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiRWFzdGVybiBCZXJpbmcgU2VhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkJheSBvZiBCaXNjYXkiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiRW5nbGlzaCBDaGFubmVsIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkd1bGYgb2YgTWV4aWNvIFN1bW1lciIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJHdWxmIG9mIEFsYXNrYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJHcmVlbmxhbmQiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTiBHdWxmIG9mIFN0LiBMYXdyZW5jZSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJTIEd1bGYgb2YgU3QuIExhd3JlbmNlIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkljZWxhbmQiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiSXJpc2ggU2VhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk1lZGl0ZXJyYW5lYW4iLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTmFtaWJpYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJORSBVUyBGYWxsIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5FIFVTIFNwcmluZyIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJOIElyZWxhbmQgUTEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTiBJcmVsYW5kIFE0IiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkJhcmVudHMgU2VhIE5vcndheSBRMyIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJOIFNlYSBRMSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJOIFNlYSBRMyIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJDaGF0aGFtIFJpc2UgTloiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiRSBDb2FzdCBTIElzbGFuZCBOWiIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJXIENvYXN0IFMgSXNsYW5kIE5aIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlBvcnR1Z2FsIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlMgR2VvcmdpYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiU2NvdGlhbiBTaGVsZiBTdW1tZXIiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlNFIFVTIEZhbGwiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlNFIFVTIFNwcmluZyIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiU0UgVVMgU3VtbWVyIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJXIENvYXN0IFVTIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJBdGxhbnRpYyBPY2VhbiBaQSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiSW5kaWFuIE9jZWFuIFpBIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiUm9ja2FsbCBQbGF0ZWF1IiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJTY290bGFuZCBTaGVsZiBTZWEgUTEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlNjb3RsYW5kIFNoZWxmIFNlYSBRNCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiRmFsa2xhbmQgSXNsYW5kcyIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiR3VsZiBvZiBNZXhpY28gRmFsbCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiQmFyZW50cyBTZWEgTm9yd2F5IFExIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJTdWItQXJjdGljIE5aIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJTY290aWFuIFNoZWxmIFNwcmluZyIpLAogICAgICAgICAgICAgICAgICAgICAgICAgIHN1cnZleV91bml0ID0gYygKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJBSSIsICAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJCSVRTLTEiLCAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJCSVRTLTQiLCAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJDSEwiLCAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJERk8tTkYiLCAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJERk8tUUNTIiwgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJFQlMiLCAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJFVkhPRSIsICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJGUi1DR0ZTIiwgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJHTUVYLVN1bW1lciIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiR09BIiwgICAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiR1JMLURFIiwgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiR1NMLU4iLCAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiR1NMLVMiLCAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiSUNFLUdGUyIsICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiSUUtSUdGUyIsICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTUVESVRTIiwgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTkFNIiwgICAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTkVVUy1GYWxsIiwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTkVVUy1TcHJpbmciLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5JR0ZTLTEiLCAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5JR0ZTLTQiLCAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5vci1CVFMtMyIsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5TLUlCVFMtMSIsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5TLUlCVFMtMyIsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5aLUNIQVQiLCAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5aLUVDU0kiLCAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5aLVdDU0kiLCAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlBULUlCVFMiLCAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlMtR0VPUkciLCAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlNDUy1TVU1NRVIiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlNFVVMtZmFsbCIsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlNFVVMtc3ByaW5nIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJTRVVTLXN1bW1lciIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiV0NBTk4iLCAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiWkFGLUFUTCIsICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiWkFGLUlORCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiUk9DS0FMTCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiU1dDLUlCVFMtMSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiU1dDLUlCVFMtNCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiRkFMSyIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiR01FWC1GYWxsIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJOb3ItQlRTLTEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5aLVNVQkEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlNDUy1TUFJJTkciCiAgICAgICAgICAgICAgICAgICAgICAgICAgKSkKCgpjb2xvcl9saW5rIDwtIG5hbWVfaGVscGVyW2NvbG9yX2xpbmssIG9uID0gInN1cnZleV91bml0Il0KCmBgYAoKCiMjTWFrZSBHQU1zCgpCcmF5IEN1cnRpcwpgYGB7ciBicmF5IGN1cnRpcyBnYW1zfQpicmF5X2N1cnRpc190b3RhbF9nYW0gPC0gZ2FtKGJyYXlfY3VydGlzX2Rpc3NpbWlsYXJpdHlfdG90YWxfbWVhbiB+IHllYXIgKyBzKHN1cnZleV91bml0LCB5ZWFyLCBicyA9ICJmcyIsIG0gPSAxKSwjcmFuZG9tIHNtb290aAogICAgICAgICAgICAgICAgICAgICAgICAgICAgZGF0YSA9IGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMucikKCmBgYAoKCiMjTWFrZSBMTUVSUwoKQnJheQoqVGhlc2UgYWxsIGNvbnZlcmdlZCoKYGBge3IgYnJheX0KI3J1bm5pbmcgd2l0aCBsbWUgaW5zdGVhZCBvZiBsbWVyIGdhdmUgc2FtZSByZXN1bHRzLCBidXQgYWxsb3dlZCBmb3IgY2FsY3VsYXRpb24gb2YgcC12YWx1ZQpicmF5X2N1cnRpc190b3RhbF9sbWUgPC0gbG1lKGJyYXlfY3VydGlzX2Rpc3NpbWlsYXJpdHlfdG90YWxfbWVhbiB+IHllYXJfYWRqLCByYW5kb20gPSAofjEgKyB5ZWFyX2FkanxzdXJ2ZXlfdW5pdCksZGF0YSA9IGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMucikKCiNidXQgYWxzbyBydW4gd2l0aCBsbWVyIGZvciBjb25maW50CiN0b3RhbApicmF5X2N1cnRpc190b3RhbF9sbWVyIDwtIGxtZXIoYnJheV9jdXJ0aXNfZGlzc2ltaWxhcml0eV90b3RhbF9tZWFuIH4geWVhcl9hZGogKyAoMSArIHllYXJfYWRqfHN1cnZleV91bml0KSxkYXRhID0gZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yKQojYmFsYW5jZWQgKGZvciBjb21wYXJpc29uIHBsb3QpCmJyYXlfY3VydGlzX2JhbGFuY2VkX2xtZXIgPC0gbG1lcihicmF5X2N1cnRpc19kaXNzaW1pbGFyaXR5X2JhbGFuY2VkX21lYW4gfiB5ZWFyX2FkaiArICgxICsgeWVhcl9hZGp8c3VydmV5X3VuaXQpLGRhdGEgPSBkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnIpCiNncmFkaWVudCAoZm9yIGNvbXBhcmlzb24gcGxvdCkKYnJheV9jdXJ0aXNfZ3JhZGllbnRfbG1lciA8LSBsbWVyKGJyYXlfY3VydGlzX2Rpc3NpbWlsYXJpdHlfZ3JhZGllbnRfbWVhbiB+IHllYXJfYWRqICsgKDEgKyB5ZWFyX2FkanxzdXJ2ZXlfdW5pdCksZGF0YSA9IGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMucikKI0phY2NhcmQgdG90YWwgKGZvciBjb21wYXJpc29uIHBsb3QpCmphY2NhcmRfdG90YWxfbG1lciA8LSBsbWVyKGphY2NhcmRfZGlzc2ltaWxhcml0eV90b3RhbF9tZWFuIH4geWVhcl9hZGogKyAoMSArIHllYXJfYWRqfHN1cnZleV91bml0KSxkYXRhID0gZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yKQoKCnN1bW1hcnkoYnJheV9jdXJ0aXNfdG90YWxfbG1lKQphbm92YShicmF5X2N1cnRpc190b3RhbF9sbWUpCgpicmF5X2N1cnRpc190b3RhbF9jb2VmcyA8LSBkYXRhLnRhYmxlKGNvZWYoYnJheV9jdXJ0aXNfdG90YWxfbG1lKSkKYnJheV9jdXJ0aXNfdG90YWxfY29lZnNbLHN1cnZleV91bml0IDo9IHJvd25hbWVzKGNvZWYoYnJheV9jdXJ0aXNfdG90YWxfbG1lKSldWyxZZWFyIDo9IHJvdW5kKHllYXJfYWRqLDUpXVssSW50ZXJjZXB0IDo9IHJvdW5kKGAoSW50ZXJjZXB0KWAsMildCiNWaWV3KGJyYXlfY3VydGlzX3RvdGFsX2NvZWZzKQoKYnJheV9jdXJ0aXNfdG90YWxfY29lZnMgPC0gYnJheV9jdXJ0aXNfdG90YWxfY29lZnNbY29sb3JfbGluaywgb24gPSAic3VydmV5X3VuaXQiXQoKYnJheV9jdXJ0aXNfdG90YWxfY29lZnMuZXhwIDwtIGJyYXlfY3VydGlzX3RvdGFsX2NvZWZzWywuKFN1cnZleV9OYW1lX1NlYXNvbiwgSW50ZXJjZXB0LCBZZWFyKV0KCiNleHBvcnQgdGhpcyB0YWJsZQpmd3JpdGUoYnJheV9jdXJ0aXNfdG90YWxfY29lZnMuZXhwLCBmaWxlID0gaGVyZTo6aGVyZSgib3V0cHV0IiwiYnJheV9jdXJ0aXNfdG90YWxfY29lZnMuZXhwLmNzdiIpKQpgYGAKCkdldCBMTUVSIG1vZGVsIGFzIHByZWRpY3Rpb25zCmBgYHtyIFRPVEFMIEJDfQoKIyBuZWVkIHRvIHNvcnQgb3V0IHllYXIgaW4gc2VxIHZlcnN1cyBvdmVyYWxsIHllYXIgbW9kZWxzCiNuZXcgZGF0YSBmb3IgbG1lcgpsbWVyX3llYXIgPC0gc2VxKG1pbihkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLHllYXJdKSwgbWF4KGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclsseWVhcl0pLCBieSA9IDAuMSkKCmxtZXJfeWVhcl9hZGogPC0gc2VxKG1pbihkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLHllYXJfYWRqXSksIG1heChkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLHllYXJfYWRqXSksIGJ5ID0gMC4xKQoKI3ByZWRpY3QgYXZlcmFnZSBsbWVyCmxtZXJfYnJheV90b3RhbF9wcmVkaWN0aW9ucyA8LSBkYXRhLnRhYmxlKHllYXIgPSBsbWVyX3llYXIsIHllYXJfYWRqID0gbG1lcl95ZWFyX2FkaikKCiNjb25maWRlbmNlIGludGVydmFscwpicmF5X2N1cnRpc190b3RhbF9sbWVyX2NvbmZpbnQgPC0gY29uZmludChicmF5X2N1cnRpc190b3RhbF9sbWVyKQoKI3BvcHVsYXRlIGRhdGEgdGFibGUgb2YgbG1lciBwcmVkaWN0aW9ucwpsbWVyX2JyYXlfdG90YWxfcHJlZGljdGlvbnNbLGJyYXlfY3VydGlzX2xtZXJfcHJlZHMgOj0gZml4ZWYoYnJheV9jdXJ0aXNfdG90YWxfbG1lcilbWzFdXSArIGZpeGVmKGJyYXlfY3VydGlzX3RvdGFsX2xtZXIpW1syXV0gKiB5ZWFyX2Fkal1bLGJyYXlfY3VydGlzX2xtZXJfcHJlZHNfbG93ZXJDSSA6PSBicmF5X2N1cnRpc190b3RhbF9sbWVyX2NvbmZpbnRbNV0gKyBicmF5X2N1cnRpc190b3RhbF9sbWVyX2NvbmZpbnRbNl0gKiB5ZWFyX2Fkal1bLGJyYXlfY3VydGlzX2xtZXJfcHJlZHNfdXBwZXJDSSA6PSBicmF5X2N1cnRpc190b3RhbF9sbWVyX2NvbmZpbnRbMTFdICsgYnJheV9jdXJ0aXNfdG90YWxfbG1lcl9jb25maW50WzEyXSAqIHllYXJfYWRqXQpgYGAKCmBgYHtyIEJBTEFOQ0VEIEJDfQoKIyBuZWVkIHRvIHNvcnQgb3V0IHllYXIgaW4gc2VxIHZlcnN1cyBvdmVyYWxsIHllYXIgbW9kZWxzCiNuZXcgZGF0YSBmb3IgbG1lcgoKCiNwcmVkaWN0IGF2ZXJhZ2UgbG1lcgpsbWVyX2JyYXlfYmFsYW5jZWRfcHJlZGljdGlvbnMgPC0gZGF0YS50YWJsZSh5ZWFyID0gbG1lcl95ZWFyLCB5ZWFyX2FkaiA9IGxtZXJfeWVhcl9hZGopCgojY29uZmlkZW5jZSBpbnRlcnZhbHMKYnJheV9jdXJ0aXNfYmFsYW5jZWRfbG1lcl9jb25maW50IDwtIGNvbmZpbnQoYnJheV9jdXJ0aXNfYmFsYW5jZWRfbG1lcikKCiNwb3B1bGF0ZSBkYXRhIHRhYmxlIG9mIGxtZXIgcHJlZGljdGlvbnMKbG1lcl9icmF5X2JhbGFuY2VkX3ByZWRpY3Rpb25zWyxicmF5X2N1cnRpc19sbWVyX3ByZWRzIDo9IGZpeGVmKGJyYXlfY3VydGlzX2JhbGFuY2VkX2xtZXIpW1sxXV0gKyBmaXhlZihicmF5X2N1cnRpc19iYWxhbmNlZF9sbWVyKVtbMl1dICogeWVhcl9hZGpdWyxicmF5X2N1cnRpc19sbWVyX3ByZWRzX2xvd2VyQ0kgOj0gYnJheV9jdXJ0aXNfYmFsYW5jZWRfbG1lcl9jb25maW50WzVdICsgYnJheV9jdXJ0aXNfYmFsYW5jZWRfbG1lcl9jb25maW50WzZdICogeWVhcl9hZGpdWyxicmF5X2N1cnRpc19sbWVyX3ByZWRzX3VwcGVyQ0kgOj0gYnJheV9jdXJ0aXNfYmFsYW5jZWRfbG1lcl9jb25maW50WzExXSArIGJyYXlfY3VydGlzX2JhbGFuY2VkX2xtZXJfY29uZmludFsxMl0gKiB5ZWFyX2Fkal0KYGBgCgpgYGB7ciBHUkFESUVOVCBCQ30KCiMgbmVlZCB0byBzb3J0IG91dCB5ZWFyIGluIHNlcSB2ZXJzdXMgb3ZlcmFsbCB5ZWFyIG1vZGVscwoKCiNwcmVkaWN0IGF2ZXJhZ2UgbG1lcgpsbWVyX2JyYXlfZ3JhZGllbnRfcHJlZGljdGlvbnMgPC0gZGF0YS50YWJsZSh5ZWFyID0gbG1lcl95ZWFyLCB5ZWFyX2FkaiA9IGxtZXJfeWVhcl9hZGopCgojY29uZmlkZW5jZSBpbnRlcnZhbHMKYnJheV9jdXJ0aXNfZ3JhZGllbnRfbG1lcl9jb25maW50IDwtIGNvbmZpbnQoYnJheV9jdXJ0aXNfZ3JhZGllbnRfbG1lcikKCiNwb3B1bGF0ZSBkYXRhIHRhYmxlIG9mIGxtZXIgcHJlZGljdGlvbnMKbG1lcl9icmF5X2dyYWRpZW50X3ByZWRpY3Rpb25zWyxicmF5X2N1cnRpc19sbWVyX3ByZWRzIDo9IGZpeGVmKGJyYXlfY3VydGlzX2dyYWRpZW50X2xtZXIpW1sxXV0gKyBmaXhlZihicmF5X2N1cnRpc19ncmFkaWVudF9sbWVyKVtbMl1dICogeWVhcl9hZGpdWyxicmF5X2N1cnRpc19sbWVyX3ByZWRzX2xvd2VyQ0kgOj0gYnJheV9jdXJ0aXNfZ3JhZGllbnRfbG1lcl9jb25maW50WzVdICsgYnJheV9jdXJ0aXNfZ3JhZGllbnRfbG1lcl9jb25maW50WzZdICogeWVhcl9hZGpdWyxicmF5X2N1cnRpc19sbWVyX3ByZWRzX3VwcGVyQ0kgOj0gYnJheV9jdXJ0aXNfZ3JhZGllbnRfbG1lcl9jb25maW50WzExXSArIGJyYXlfY3VydGlzX2dyYWRpZW50X2xtZXJfY29uZmludFsxMl0gKiB5ZWFyX2Fkal0KYGBgCgpgYGB7ciB0b3RhbCBqYWNjYXJkfQoKIyBuZWVkIHRvIHNvcnQgb3V0IHllYXIgaW4gc2VxIHZlcnN1cyBvdmVyYWxsIHllYXIgbW9kZWxzCgoKI3ByZWRpY3QgYXZlcmFnZSBsbWVyCmxtZXJfamFjY2FyZF90b3RhbF9wcmVkaWN0aW9ucyA8LSBkYXRhLnRhYmxlKHllYXIgPSBsbWVyX3llYXIsIHllYXJfYWRqID0gbG1lcl95ZWFyX2FkaikKCiNjb25maWRlbmNlIGludGVydmFscwpqYWNjYXJkX3RvdGFsX2xtZXJfY29uZmludCA8LSBjb25maW50KGphY2NhcmRfdG90YWxfbG1lcikKCiNwb3B1bGF0ZSBkYXRhIHRhYmxlIG9mIGxtZXIgcHJlZGljdGlvbnMKbG1lcl9qYWNjYXJkX3RvdGFsX3ByZWRpY3Rpb25zWyxqYWNjYXJkX2xtZXJfcHJlZHMgOj0gZml4ZWYoamFjY2FyZF90b3RhbF9sbWVyKVtbMV1dICsgZml4ZWYoamFjY2FyZF90b3RhbF9sbWVyKVtbMl1dICogeWVhcl9hZGpdWyxqYWNjYXJkX2xtZXJfcHJlZHNfbG93ZXJDSSA6PSBqYWNjYXJkX3RvdGFsX2xtZXJfY29uZmludFs1XSArIGphY2NhcmRfdG90YWxfbG1lcl9jb25maW50WzZdICogeWVhcl9hZGpdWyxqYWNjYXJkX2xtZXJfcHJlZHNfdXBwZXJDSSA6PSBqYWNjYXJkX3RvdGFsX2xtZXJfY29uZmludFsxMV0gKyBqYWNjYXJkX3RvdGFsX2xtZXJfY29uZmludFsxMl0gKiB5ZWFyX2Fkal0KYGBgCiMjI01vdmUgZm9yd2FyZCB3aXRoIEJyYXkgQ3VydGlzIHRvdGFsIChvdGhlcnMgdG8gc3VwcGxlbWVudCwgYWx0aG91Z2ggd2lsbCBpbmNsdWRlIGluIHBsb3QgYmVsb3cpCgoKQ29lZmZpY2llbnRzIGZvciBMTUVSIGJ5IHN1cnZleV91bml0CmBgYHtyfQojdW5pcXVlIHN1cnZleSB1bml0IHllYXIgY29tYm9zCnN1cnZleV91bml0X3NhbXBsaW5nX3llYXJzIDwtIHVuaXF1ZShkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLC4oc3VydmV5X3VuaXQsIHllYXJfYWRqLCB5ZWFyLCB5ZWFyc19zYW1wbGVkKV0pCgojIHNlZSBncm91cCBjb2VmZmljaWVudHMKbW9kZWxfY29lZnNfcmVkdWNlZCA8LSBkYXRhLnRhYmxlKHRyYW5zZm9ybShhcy5kYXRhLmZyYW1lKHJhbmVmKGJyYXlfY3VydGlzX3RvdGFsX2xtZXIpKSwgbHdyID0gY29uZHZhbCAtIDEuOTYqY29uZHNkLCB1cHIgPSBjb25kdmFsICsgMS45Nipjb25kc2QpKQojaHR0cHM6Ly9zdGFja292ZXJmbG93LmNvbS9xdWVzdGlvbnMvNjk4MDU1MzIvZXh0cmFjdC10aGUtY29uZmlkZW5jZS1pbnRlcnZhbHMtb2YtbG1lci1yYW5kb20tZWZmZWN0cy1wbG90dGVkLXdpdGgtZG90cGxvdHJhCgoKI09OTFkgU0xPUEVTCm1vZGVsX2NvZWZzX3JlZHVjZWQgPC0gbW9kZWxfY29lZnNfcmVkdWNlZFt0ZXJtID09ICJ5ZWFyX2FkaiIsXQoKbW9kZWxfY29lZnNfcmVkdWNlZFssc3VydmV5X3VuaXQgOj0gZ3JwXVsseWVhcl9hZGogOj0gY29uZHZhbF0KCiNtZXJnZSB3aXRoIGR1cmF0aW9uIG9mIHN1cnZleQptb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aCA8LSBtb2RlbF9jb2Vmc19yZWR1Y2VkW3N1cnZleV91bml0X3NhbXBsaW5nX3llYXJzLCBvbiA9ICJzdXJ2ZXlfdW5pdCJdCgoKbW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGhbLHllYXJzX3NhbXBsZWQgOj0gYXMubnVtZXJpYyh5ZWFyc19zYW1wbGVkKV1bLERpcmVjdGlvbmFsX0NoYW5nZSA6PSBpZmVsc2UoeWVhcl9hZGogPiAwLCAiRGlmZmVyZW50aWF0aW9uIiwiSG9tb2dlbml6YXRpb24iKV0KCiNkb2VzIGl0IGNyb3NzIHplcm8/Cm1vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoWyxzaWduaWZpY2FudCA6PSBpZmVsc2UobHdyID4wICYgdXByPjAsVCxpZmVsc2UobHdyPDAgJiB1cHI8MCxULEYpKV0KCiNzaWduaWZpY2FudCBkaXJlY3Rpb25hbCBjaGFuZ2UKbW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGhbLERpcmVjdGlvbmFsX0NoYW5nZV9zaWcgOj0gaWZlbHNlKHllYXJfYWRqID4gMCAmIHNpZ25pZmljYW50ID09IFQsICJEaWZmZXJlbnRpYXRpb24iLGlmZWxzZSh5ZWFyX2FkaiA8IDAgJiBzaWduaWZpY2FudCA9PSBULCAiSG9tb2dlbml6YXRpb24iLCAiTm8gdHJlbmQgaW5cbmRpc3NpbWlsYXJpdHkiKSldCgoKI21pbiBtYXggZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllcwptaW5fYnJheV9yZWR1Y2VkIDwtIG1pbihkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLGJyYXlfY3VydGlzX2Rpc3NpbWlsYXJpdHlfdG90YWxfbWVhbl0sIG5hLnJtID0gVCkKbWF4X2JyYXlfcmVkdWNlZCA8LSBtYXgoZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyxicmF5X2N1cnRpc19kaXNzaW1pbGFyaXR5X3RvdGFsX21lYW5dLCBuYS5ybSA9IFQpCgptb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aCA8LSBtb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aFtjb2xvcl9saW5rLCBvbiA9ICJzdXJ2ZXlfdW5pdCJdCgojZGVsZXRlIGFueSBOQXMKbW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGggPC0gbmEub21pdChtb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aCwgY29scyA9ICJzaWduaWZpY2FudCIpCgojb3JkZXIgdGFibGUgYnkgY29lZmZpY2llbnQKc2V0b3JkZXIobW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGgsIHllYXJfYWRqKQoKQkNfdG90YWxfbW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGgudW5pcXVlIDwtIHVuaXF1ZShtb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aFssLihjb25kdmFsLGNvbmRzZCwgbHdyLCB1cHIsIHN1cnZleV91bml0LCB5ZWFyX2FkaiwgeWVhcnNfc2FtcGxlZCwgRGlyZWN0aW9uYWxfQ2hhbmdlLCBoZXgsIFN1cnZleV9OYW1lX1NlYXNvbiwgc2lnbmlmaWNhbnQsIERpcmVjdGlvbmFsX0NoYW5nZV9zaWcpXSkgCgojZXh0cmFjdCBjb2xvciBoZXhlcwojeWVhciBhZGogY29lZiBvcmRlcgpjb2xvcl95ZWFyX2Fkal9vcmRlciA8LSBCQ190b3RhbF9tb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aC51bmlxdWVbLGhleF0KCiNhbHBoYWJldGljYWwgb3JkZXIKQkNfdG90YWxfbW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGgudW5pcXVlLmFscGhhIDwtIHNldG9yZGVyKEJDX3RvdGFsX21vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoLnVuaXF1ZSwgU3VydmV5X05hbWVfU2Vhc29uKQoKQkNfdG90YWxfbW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGgudW5pcXVlLmFscGhhWyx0cmVuZF9jb2xvciA6PSBpZmVsc2UoRGlyZWN0aW9uYWxfQ2hhbmdlX3NpZyA9PSAiSG9tb2dlbml6YXRpb24iLCAiI2U3YWM1YiIsIGlmZWxzZShEaXJlY3Rpb25hbF9DaGFuZ2Vfc2lnID09ICJEaWZmZXJlbnRpYXRpb24iLCIjOTFjODc0IiwiI2NiYmZkZSIpKV0KCmNvbG9yX2FscGhhX29yZGVyIDwtIEJDX3RvdGFsX21vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoLnVuaXF1ZS5hbHBoYVssaGV4XQpjb2xvcl9hbHBoYV9vcmRlcl9ieXRyZW5kIDwtIEJDX3RvdGFsX21vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoLnVuaXF1ZS5hbHBoYVssIHRyZW5kX2NvbG9yXQoKc2F2ZVJEUyhCQ190b3RhbF9tb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aC51bmlxdWUsIGhlcmU6OmhlcmUoIm91dHB1dCIsInJlZ2lvbl9zdGF0cyIsIkJDX3RvdGFsX21vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoLnVuaXF1ZS5SZHMiKSkKYGBgCgpUb3RhbCB2ZXJzdXMgYmFsYW5jZWQgQkMgcGxvdApgYGB7cn0KZ2dwbG90KGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMucikgKwogIGdlb21fcG9pbnQoYWVzKGJyYXlfY3VydGlzX2Rpc3NpbWlsYXJpdHlfdG90YWxfbWVhbiwgYnJheV9jdXJ0aXNfZGlzc2ltaWxhcml0eV9iYWxhbmNlZF9tZWFuKSkgKwogIGdlb21fYWJsaW5lKHNsb3BlID0gMSwgaW50ZXJjcGV0ID0gMCkgKwogIGdlb21fc21vb3RoKGFlcyhicmF5X2N1cnRpc19kaXNzaW1pbGFyaXR5X3RvdGFsX21lYW4sIGJyYXlfY3VydGlzX2Rpc3NpbWlsYXJpdHlfYmFsYW5jZWRfbWVhbikpICsKICB0aGVtZV9jbGFzc2ljKCkgKwogIGxhYnMoeCA9ICJUb3RhbCBCQyBkaXNzaW1pbGFyaXR5IiwgIHkgPSAiQmFsYW5jZWQgY2hhbmdlcyBpbiBhYnVuZGFuY2UvYmlvbWFzcyIpCgpnZ3Bsb3QoZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yKSArCiAgZ2VvbV9wb2ludChhZXMoYnJheV9jdXJ0aXNfZGlzc2ltaWxhcml0eV90b3RhbF9tZWFuLCBicmF5X2N1cnRpc19kaXNzaW1pbGFyaXR5X2dyYWRpZW50X21lYW4pKSArCiAgZ2VvbV9hYmxpbmUoc2xvcGUgPSAxLCBpbnRlcmNwZXQgPSAwKSArCiAgZ2VvbV9zbW9vdGgoYWVzKGJyYXlfY3VydGlzX2Rpc3NpbWlsYXJpdHlfdG90YWxfbWVhbiwgYnJheV9jdXJ0aXNfZGlzc2ltaWxhcml0eV9ncmFkaWVudF9tZWFuKSkgKwogIHRoZW1lX2NsYXNzaWMoKSArCiAgbGFicyh4ID0gIlRvdGFsIEJDIGRpc3NpbWlsYXJpdHkiLCAgeSA9ICJBYnVuZGFuY2UvYmlvbWFzcyBncmFkaWVudCIpCgojYWxsIG9uIG9uZSBwbG90CiNyZWR1Y2VkCmRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuci5yIDwtIGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclssLih5ZWFyLCBicmF5X2N1cnRpc19kaXNzaW1pbGFyaXR5X3RvdGFsX21lYW4sIGJyYXlfY3VydGlzX2Rpc3NpbWlsYXJpdHlfZ3JhZGllbnRfbWVhbiwgYnJheV9jdXJ0aXNfZGlzc2ltaWxhcml0eV9iYWxhbmNlZF9tZWFuLCBqYWNjYXJkX2Rpc3NpbWlsYXJpdHlfdG90YWxfbWVhbildCgojd2lkZSB0byBsb25nCmRpc3NpbV9sb25nIDwtIG1lbHQoZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yLnIsIGlkLnZhcnMgPSBjKCJ5ZWFyIiksIHZhcmlhYmxlLm5hbWUgPSAiRGlzc2ltaWxhcml0eSBtZXRyaWMiKQoKZGlzc2ltX2xvbmdbLGBEaXNzaW1pbGFyaXR5IG1ldHJpY2AgOj0gZmFjdG9yKGBEaXNzaW1pbGFyaXR5IG1ldHJpY2AsIGxldmVscyA9IGMoICJicmF5X2N1cnRpc19kaXNzaW1pbGFyaXR5X3RvdGFsX21lYW4iLCJicmF5X2N1cnRpc19kaXNzaW1pbGFyaXR5X2JhbGFuY2VkX21lYW4iLCJicmF5X2N1cnRpc19kaXNzaW1pbGFyaXR5X2dyYWRpZW50X21lYW4iLCAiamFjY2FyZF9kaXNzaW1pbGFyaXR5X3RvdGFsX21lYW4iKSldCgpCQ19kaXNzaW1pbGFyaXR5X2NvbXBvbmVudHMgPC0gZ2dwbG90KCkgKwogIGdlb21fcG9pbnQoZGF0YSA9IGRpc3NpbV9sb25nW2BEaXNzaW1pbGFyaXR5IG1ldHJpY2AgIT0gImphY2NhcmRfZGlzc2ltaWxhcml0eV90b3RhbF9tZWFuIixdLCBhZXMoeCA9IHllYXIsIHkgPSB2YWx1ZSwgY29sb3IgPSBgRGlzc2ltaWxhcml0eSBtZXRyaWNgKSwgYWxwaGEgPSAwLjIsIHNpemUgPSAxKSArCiAgZ2VvbV9yaWJib24oZGF0YSA9IGxtZXJfYnJheV90b3RhbF9wcmVkaWN0aW9ucywgYWVzKHggPSB5ZWFyLCB5bWluID0gYnJheV9jdXJ0aXNfbG1lcl9wcmVkc19sb3dlckNJLCB5bWF4ID0gYnJheV9jdXJ0aXNfbG1lcl9wcmVkc191cHBlckNJKSwgZmlsbCA9ICJibGFjayIsIGFscGhhID0gMC4zKSArCiAgZ2VvbV9saW5lKGRhdGEgPSBsbWVyX2JyYXlfdG90YWxfcHJlZGljdGlvbnMsIGFlcyh4ID0geWVhciwgeSA9IGJyYXlfY3VydGlzX2xtZXJfcHJlZHMpLCBjb2xvciA9ICJibGFjayIpICsKICAjYmFsYW5jZWQKICAgIGdlb21fcmliYm9uKGRhdGEgPSBsbWVyX2JyYXlfZ3JhZGllbnRfcHJlZGljdGlvbnMsIGFlcyh4ID0geWVhciwgeW1pbiA9IGJyYXlfY3VydGlzX2xtZXJfcHJlZHNfbG93ZXJDSSwgeW1heCA9IGJyYXlfY3VydGlzX2xtZXJfcHJlZHNfdXBwZXJDSSksIGZpbGwgPSAiY2FkZXRibHVlMyIsIGFscGhhID0gMC4zKSArCiAgZ2VvbV9saW5lKGRhdGEgPSBsbWVyX2JyYXlfZ3JhZGllbnRfcHJlZGljdGlvbnMsIGFlcyh4ID0geWVhciwgeSA9IGJyYXlfY3VydGlzX2xtZXJfcHJlZHMpLCBjb2xvciA9ICJjYWRldGJsdWUzIikgKwogICNncmFkaWVudAogICAgZ2VvbV9yaWJib24oZGF0YSA9IGxtZXJfYnJheV9iYWxhbmNlZF9wcmVkaWN0aW9ucywgYWVzKHggPSB5ZWFyLCB5bWluID0gYnJheV9jdXJ0aXNfbG1lcl9wcmVkc19sb3dlckNJLCB5bWF4ID0gYnJheV9jdXJ0aXNfbG1lcl9wcmVkc191cHBlckNJKSwgZmlsbCA9ICJkYXJrc2FsbW9uIiwgYWxwaGEgPSAwLjMpICsKICBnZW9tX2xpbmUoZGF0YSA9IGxtZXJfYnJheV9iYWxhbmNlZF9wcmVkaWN0aW9ucywgYWVzKHggPSB5ZWFyLCB5ID0gYnJheV9jdXJ0aXNfbG1lcl9wcmVkcyksIGNvbG9yID0gImRhcmtzYWxtb24iKSArCiAgI2ZpeC4gY29sb3JzIG9mIHBvaW50cwogIHNjYWxlX2NvbG9yX21hbnVhbCh2YWx1ZXMgPSBjKCJibGFjayIsImRhcmtzYWxtb24iLCJjYWRldGJsdWUzIiksIGxhYmVscyA9IGMoIlRvdGFsIiwiQmFsYW5jZWQgdmFyaWF0aW9uXG5pbiBhYnVuZGFuY2UiLCJCaW9tYXNzIGdyYWRpZW50IikpICsKICBsYWJzKHggPSAiWWVhciIsIHkgPSAiVmFsdWUiKSArCiAgbGltcyh5ID0gYyhtaW4oZGlzc2ltX2xvbmckdmFsdWUpLCBtYXgoZGlzc2ltX2xvbmckdmFsdWUpKSkgKwogIHRoZW1lX2NsYXNzaWMoKSArCiAgdGhlbWUobGVnZW5kLnBvc2l0aW9uID0gYygwLjE3LDAuMzcpKQoKZ2dzYXZlKEJDX2Rpc3NpbWlsYXJpdHlfY29tcG9uZW50cywgcGF0aCA9IGhlcmU6OmhlcmUoImZpZ3VyZXMiKSwgZmlsZW5hbWU9KCJCQ19kaXNzaW1pbGFyaXR5X2NvbXBvbmVudHMuanBnIiksIHdpZHRoID0gNi41LCBoZWlnaHQgPSA1LCB1bml0cyA9ICJpbiIpCgojSmFjY2FyZApKYWNjYXJkX2Rpc3NpbWlsYXJpdHkgPC0gZ2dwbG90KCkgKwogIGdlb21fcG9pbnQoZGF0YSA9IGRpc3NpbV9sb25nW2BEaXNzaW1pbGFyaXR5IG1ldHJpY2AgPT0gImphY2NhcmRfZGlzc2ltaWxhcml0eV90b3RhbF9tZWFuIl0sIGFlcyh4ID0geWVhciwgeSA9IHZhbHVlKSwgYWxwaGEgPSAwLjIsIHNpemUgPSAxKSArCiAgZ2VvbV9yaWJib24oZGF0YSA9IGxtZXJfamFjY2FyZF90b3RhbF9wcmVkaWN0aW9ucywgYWVzKHggPSB5ZWFyLCB5bWluID0gamFjY2FyZF9sbWVyX3ByZWRzX2xvd2VyQ0ksIHltYXggPSBqYWNjYXJkX2xtZXJfcHJlZHNfdXBwZXJDSSksIGZpbGwgPSAiYmxhY2siLCBhbHBoYSA9IDAuMykgKwogIGdlb21fbGluZShkYXRhID0gbG1lcl9qYWNjYXJkX3RvdGFsX3ByZWRpY3Rpb25zLCBhZXMoeCA9IHllYXIsIHkgPSBqYWNjYXJkX2xtZXJfcHJlZHMpLCBjb2xvciA9ICJibGFjayIpICsKICBsYWJzKHggPSAiWWVhciIsIHkgPSAiVG90YWwgSmFjY2FyZCBkaXNzaW1pbGFyaXR5IikgKwogIGxpbXMoeSA9IGMobWluKGRpc3NpbV9sb25nJHZhbHVlKSwgbWF4KGRpc3NpbV9sb25nJHZhbHVlKSkpICsKICB0aGVtZV9jbGFzc2ljKCkKCmdnc2F2ZShKYWNjYXJkX2Rpc3NpbWlsYXJpdHksIHBhdGggPSBoZXJlOjpoZXJlKCJmaWd1cmVzIiksIGZpbGVuYW1lPSgiSmFjY2FyZF9kaXNzaW1pbGFyaXR5LmpwZyIpLCB3aWR0aCA9IDYuNSwgaGVpZ2h0ID0gNSwgdW5pdHMgPSAiaW4iKQoKI21lcmdlCmxpYnJhcnkoY293cGxvdCkKYWxsX2Rpc3NpbWlsYXJpdGllc19vdmVyX3RpbWUgPC0gcGxvdF9ncmlkKEJDX2Rpc3NpbWlsYXJpdHlfY29tcG9uZW50cyArIHRoZW1lKGxlZ2VuZC5iYWNrZ3JvdW5kID0gZWxlbWVudF9yZWN0KGZpbGwgPSAidHJhbnNwYXJlbnQiKSksCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBKYWNjYXJkX2Rpc3NpbWlsYXJpdHksCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBuY29sID0gMiwgYWxpZ24gPSAiaHYiLCBsYWJlbHMgPSBjKCJhLiIsImIuIikpCgpnZ3NhdmUoYWxsX2Rpc3NpbWlsYXJpdGllc19vdmVyX3RpbWUsIGZpbGVuYW1lID0gImFsbF9kaXNzaW1pbGFyaXRpZXNfb3Zlcl90aW1lLmpwZyIsIHBhdGggPSBoZXJlOjpoZXJlKCJmaWd1cmVzIiksIHdpZHRoID0gMTAsIGhlaWdodCA9IDQuNSwgdW5pdHMgPSAiaW4iKQoKICAKYGBgCgoKCkJhciBQbG90IENvZWZmaWNpZW50IExNRVIKYGBge3IgYmFyIHBsb3Qgb2YgY29lZmZpY2llbnRzfQpCQ190b3RhbF9EaXNzaW1pbGFyaXR5X0NvZWZfZXJyb3JiYXJfcmVkdWNlZCA8LSBnZ3Bsb3QoKSArCiAgICBnZW9tX2Vycm9yYmFyKGRhdGEgPSBtb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aCwgYWVzKHggPSByZW9yZGVyKFN1cnZleV9OYW1lX1NlYXNvbiwgeWVhcl9hZGopICwgeSA9IHllYXJfYWRqLCBsYWJlbCA9IFN1cnZleV9OYW1lX1NlYXNvbiwgeW1pbiA9IGx3ciwgeW1heCA9IHVwciksIGZpbGwgPSAiZ3JleSIsIHdpZHRoID0gMCkgKyAjYWRkIGNvbmZpZGVuY2UgaW50ZXJ2YWxzCiAgZ2VvbV9wb2ludChkYXRhID0gbW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGgsIGFlcyh4ID0gcmVvcmRlcihTdXJ2ZXlfTmFtZV9TZWFzb24sIHllYXJfYWRqKSAsIHkgPSB5ZWFyX2FkaiwgbGFiZWwgPSBTdXJ2ZXlfTmFtZV9TZWFzb24sIHNpemUgPSB5ZWFyc19zYW1wbGVkLCBmaWxsID0gRGlyZWN0aW9uYWxfQ2hhbmdlX3NpZywgY29sb3IgPSBEaXJlY3Rpb25hbF9DaGFuZ2Vfc2lnKSwgc3RhdCA9ICdpZGVudGl0eScsIHNoYXBlID0gMjEpICsKICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSBjKCJ3aGl0ZSIsImJsYWNrIiwiZ3JleSIpLCBuYW1lID0gIkRpc3NpbWlsYXJpdHkgdHJlbmQiLCBndWlkZT0ibm9uZSIpICsKICBzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzID0gYygiYmxhY2siLCJibGFjayIsImdyZXkiKSwgbmFtZSA9ICJEaXNzaW1pbGFyaXR5IHRyZW5kIiwgZ3VpZGU9Im5vbmUiKSArCiAgc2NhbGVfc2l6ZV9iaW5uZWQocmFuZ2UgPSBjKDEsOCksIG5hbWUgPSAiU3VydmV5IHBlcmlvZCBsZW5ndGhcbih5ZWFycykiKSArCiAgZ2VvbV9obGluZSh5aW50ZXJjZXB0ID0gMCkgKwogIHNjYWxlX3lfY29udGludW91cyhicmVha3MgPSBzZXEoLTAuMDA1LCAwLjAwNzUsIGJ5ID0gMC4wMDI1KSwgbGFiZWxzID0gYygiLTAuMDA1IiwiIiwiMCIsICIiLCAiMC4wMDUiLCAgIiIpKSArCiAgeGxhYigiU3VydmV5IHVuaXQiKSArCiAgeWxhYigizrItZGl2ZXJzaXR5IHRyZW5kIikgKyAjdG90YWwgQkMgZGlzc2ltaWxhcml0eSB0cmVuZAogIGNvb3JkX2ZsaXAoKSArCiAgdGhlbWVfY2xhc3NpYygpICsKICB0aGVtZShheGlzLnRleHQueSA9IGVsZW1lbnRfdGV4dChjb2xvdXIgPSBjb2xvcl95ZWFyX2Fkal9vcmRlciwgZmFjZSA9ICJib2xkIiksIGF4aXMudGl0bGUueSA9IGVsZW1lbnRfYmxhbmsoKSwgYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSwgYXhpcy50aXRsZS54ID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNSksIGxlZ2VuZC5wb3NpdGlvbiA9IGMoMC4yNSwwLjcpLCBsZWdlbmQuZGlyZWN0aW9uID0gInZlcnRpY2FsIikKCiNwdWxsIGxlZ2VuZCBmb3IgaG9tb2dlbml6YXRpb24KZGlyZWN0aW9uYWxfY2hhbmdlX2xlZ2VuZF9wbG90IDwtIEJDX3RvdGFsX0Rpc3NpbWlsYXJpdHlfQ29lZl9lcnJvcmJhcl9yZWR1Y2VkICsgCiAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gYygid2hpdGUiLCJibGFjayIsImdyZXkiKSwgbmFtZSA9ICJEaXNzaW1pbGFyaXR5IHRyZW5kIikgKwogIHNjYWxlX2NvbG9yX21hbnVhbCh2YWx1ZXMgPSBjKCJibGFjayIsImJsYWNrIiwiZ3JleSIpLCBuYW1lID0gIkRpc3NpbWlsYXJpdHkgdHJlbmQiKSArCiAgc2NhbGVfc2l6ZV9iaW5uZWQocmFuZ2UgPSBjKDEsOCksIG5hbWUgPSAiU3VydmV5IHBlcmlvZCBsZW5ndGgiLCBndWlkZSA9ICJub25lIikgKwogIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbiA9ICJyaWdodCIsIGxlZ2VuZC5iYWNrZ3JvdW5kID0gZWxlbWVudF9yZWN0KGZpbGw9ICJ0cmFuc3BhcmVudCIpKSArCiAgZ3VpZGVzKGNvbG91ciA9IGd1aWRlX2xlZ2VuZChvdmVycmlkZS5hZXMgPSBsaXN0KHNpemU9NikpKQoKCkJDX3RvdGFsX0Rpc3NpbWlsYXJpdHlfQ29lZl9lcnJvcmJhcl9yZWR1Y2VkCgojQUxUIGdyZXkgc2NhbGUKQkNfdG90YWxfRGlzc2ltaWxhcml0eV9Db2VmX2Vycm9yYmFyX3JlZHVjZWRfZ3JleXNjYWxlIDwtIGdncGxvdCgpICsKICAgIGdlb21fZXJyb3JiYXIoZGF0YSA9IG1vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoLCBhZXMoeCA9IHJlb3JkZXIoU3VydmV5X05hbWVfU2Vhc29uLCB5ZWFyX2FkaikgLCB5ID0geWVhcl9hZGosIGxhYmVsID0gU3VydmV5X05hbWVfU2Vhc29uLCB5bWluID0gbHdyLCB5bWF4ID0gdXByKSwgZmlsbCA9ICJncmV5Iiwgd2lkdGggPSAwKSArICNhZGQgY29uZmlkZW5jZSBpbnRlcnZhbHMKICBnZW9tX3BvaW50KGRhdGEgPSBtb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aCwgYWVzKHggPSByZW9yZGVyKFN1cnZleV9OYW1lX1NlYXNvbiwgeWVhcl9hZGopICwgeSA9IHllYXJfYWRqLCBsYWJlbCA9IFN1cnZleV9OYW1lX1NlYXNvbiwgc2l6ZSA9IHllYXJzX3NhbXBsZWQsIGZpbGwgPSBEaXJlY3Rpb25hbF9DaGFuZ2Vfc2lnLCBjb2xvciA9IERpcmVjdGlvbmFsX0NoYW5nZV9zaWcpLCBzdGF0ID0gJ2lkZW50aXR5Jywgc2hhcGUgPSAyMSkgKwogIHNjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcyA9IGMoIndoaXRlIiwiYmxhY2siLCJncmV5IiksIG5hbWUgPSAiRGlzc2ltaWxhcml0eSB0cmVuZCIsIGd1aWRlPSJub25lIikgKwogIHNjYWxlX2NvbG9yX21hbnVhbCh2YWx1ZXMgPSBjKCJibGFjayIsImJsYWNrIiwiZ3JleSIpLCBuYW1lID0gIkRpc3NpbWlsYXJpdHkgdHJlbmQiLCBndWlkZT0ibm9uZSIpICsKICBzY2FsZV9zaXplX2Jpbm5lZChyYW5nZSA9IGMoMSw4KSwgbmFtZSA9ICJTdXJ2ZXkgcGVyaW9kIGxlbmd0aFxuKHllYXJzKSIpICsKICBnZW9tX2hsaW5lKHlpbnRlcmNlcHQgPSAwKSArCiAgc2NhbGVfeV9jb250aW51b3VzKGJyZWFrcyA9IHNlcSgtMC4wMDUsIDAuMDA3NSwgYnkgPSAwLjAwMjUpLCBsYWJlbHMgPSBjKCItMC4wMDUiLCIiLCIwIiwgIiIsICIwLjAwNSIsICAiIikpICsKICB4bGFiKCJTdXJ2ZXkgdW5pdCIpICsKICB5bGFiKCLOsi1kaXZlcnNpdHkgdHJlbmQiKSArICN0b3RhbCBCQyBkaXNzaW1pbGFyaXR5IHRyZW5kCiAgY29vcmRfZmxpcCgpICsKICB0aGVtZV9jbGFzc2ljKCkgKwogIHRoZW1lKGF4aXMudGV4dC55ID0gZWxlbWVudF90ZXh0KGZhY2UgPSAiYm9sZCIpLCBheGlzLnRpdGxlLnkgPSBlbGVtZW50X2JsYW5rKCksIGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNSksIGF4aXMudGl0bGUueCA9IGVsZW1lbnRfdGV4dChzaXplID0gMTUpLCBsZWdlbmQucG9zaXRpb24gPSBjKDAuMjUsMC43KSwgbGVnZW5kLmRpcmVjdGlvbiA9ICJ2ZXJ0aWNhbCIpCmBgYAoKQWx0ZXJuYXRpdmVseSwgd2UgY29sb3IgIHRoaXMgcGxvdCBieSB0cmVuZCBleHBlcmllbmNlZAoKYGBge3J9CiMiIzczQkE0RCIsIiNFMDk2MkMiLCIjY2JiZmRlIgoKQkNfdG90YWxfRGlzc2ltaWxhcml0eV9Db2VmX2Vycm9yYmFyX3JlZHVjZWRfY29sb3JieXRyZW5kIDwtIGdncGxvdCgpICsKICAgIGdlb21fZXJyb3JiYXIoZGF0YSA9IG1vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoLCBhZXMoeCA9IHJlb3JkZXIoU3VydmV5X05hbWVfU2Vhc29uLCB5ZWFyX2FkaikgLCB5ID0geWVhcl9hZGosIGxhYmVsID0gU3VydmV5X05hbWVfU2Vhc29uLCB5bWluID0gbHdyLCB5bWF4ID0gdXByKSwgZmlsbCA9ICJncmV5Iiwgd2lkdGggPSAwKSArICNhZGQgY29uZmlkZW5jZSBpbnRlcnZhbHMKICBnZW9tX3BvaW50KGRhdGEgPSBtb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aCwgYWVzKHggPSByZW9yZGVyKFN1cnZleV9OYW1lX1NlYXNvbiwgeWVhcl9hZGopICwgeSA9IHllYXJfYWRqLCBsYWJlbCA9IFN1cnZleV9OYW1lX1NlYXNvbiwgc2l6ZSA9IHllYXJzX3NhbXBsZWQsIGNvbG9yID0gRGlyZWN0aW9uYWxfQ2hhbmdlX3NpZyksIHN0YXQgPSAnaWRlbnRpdHknKSArCiAgc2NhbGVfY29sb3JfbWFudWFsKHZhbHVlcyA9IGMoIiM3M0JBNEQiLCIjRTA5NjJDIiwiI2NiYmZkZSIpLCBuYW1lID0gIkRpc3NpbWlsYXJpdHkgdHJlbmQiLCBndWlkZT0ibm9uZSIpICsKICBzY2FsZV9zaXplX2Jpbm5lZChyYW5nZSA9IGMoMSw4KSwgbmFtZSA9ICJTdXJ2ZXkgcGVyaW9kIGxlbmd0aFxuKHllYXJzKSIpICsKICBnZW9tX2hsaW5lKHlpbnRlcmNlcHQgPSAwKSArCiAgc2NhbGVfeV9jb250aW51b3VzKGJyZWFrcyA9IHNlcSgtMC4wMDUsIDAuMDA3NSwgYnkgPSAwLjAwMjUpLCBsYWJlbHMgPSBjKCItMC4wMDUiLCIiLCIwIiwgIiIsICIwLjAwNSIsICAiIikpICsKICB4bGFiKCJTdXJ2ZXkgdW5pdCIpICsKICB5bGFiKCLOsi1kaXZlcnNpdHkgdHJlbmQiKSArICN0b3RhbCBCQyBkaXNzaW1pbGFyaXR5IHRyZW5kCiAgY29vcmRfZmxpcCgpICsKICB0aGVtZV9jbGFzc2ljKCkgKwogIHRoZW1lKGF4aXMudGV4dC55ID0gZWxlbWVudF90ZXh0KGZhY2UgPSAiYm9sZCIpLCBheGlzLnRpdGxlLnkgPSBlbGVtZW50X2JsYW5rKCksIGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNSksIGF4aXMudGl0bGUueCA9IGVsZW1lbnRfdGV4dChzaXplID0gMTUpLCBsZWdlbmQucG9zaXRpb24gPSBjKDAuMywwLjgpLCBsZWdlbmQuZGlyZWN0aW9uID0gInZlcnRpY2FsIiwgbGVnZW5kLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSwgbGVnZW5kLnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNikpCgpkaXJlY3Rpb25hbF9jaGFuZ2VfbGVnZW5kX3Bsb3RfY29sb3JieXRyZW5kIDwtIEJDX3RvdGFsX0Rpc3NpbWlsYXJpdHlfQ29lZl9lcnJvcmJhcl9yZWR1Y2VkX2NvbG9yYnl0cmVuZCArIAogIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbiA9ICJyaWdodCIsIGxlZ2VuZC5iYWNrZ3JvdW5kID0gZWxlbWVudF9yZWN0KGZpbGw9ICJ0cmFuc3BhcmVudCIpLCAKICAgICAgICAgbGVnZW5kLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSwgbGVnZW5kLnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNikpICsKICBndWlkZXMoY29sb3VyID0gZ3VpZGVfbGVnZW5kKG92ZXJyaWRlLmFlcyA9IGxpc3Qoc2l6ZT02KSksIHNpemUgPSAibm9uZSIpCmBgYAoKCldhdnkgTGluZSBQbG90IGZvciBHQU1zCgpHZW5lcmF0ZSBwcmVkaWN0ZWQgdmFsdWVzCmBgYHtyIGdlbmVyYXRlIHByZWRpY3RlZCB2YWx1ZXMgR0FNfQoKI2FkZCBjb2xvcnMgYW5kIG5hbWVzIHRvIGZ1bGwgZGlzc2ltaWxhcml0eSBkYXRhIHRhYmxlCmRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuciA8LSBkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbY29sb3JfbGluaywgb24gPSAic3VydmV5X3VuaXQiXQoKI2dlbmVyYXRlIG5ldyBkYXRhIHRvIHNtb290aCBsaW5lcyAobmVlZCB5ZWFyIGFuZCBzZWFzb24gc3VydmV5IGNvbWJpbmF0aW9ucykKeWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHQgPC0gZGF0YS50YWJsZShzdXJ2ZXlfdW5pdCA9IGFzLmNoYXJhY3RlcihOVUxMKSwgeWVhciA9IGFzLm51bWVyaWMoTlVMTCksIHllYXJfYWRqID0gYXMubnVtZXJpYyhOVUxMICkpCgpmb3IgKGkgaW4gMTpsZW5ndGgoc3VydmV5X3VuaXQubGlzdCkpIHsKICAjZ2VuZXJhdGUgeWVhciB2ZWN0b3JzCiAgc3VydmV5X3VuaXRfeWVhcnMgPC0gdW5pcXVlKGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMucltzdXJ2ZXlfdW5pdCA9PSBzdXJ2ZXlfdW5pdC5saXN0W2ldLC4oc3VydmV5X3VuaXQsIHllYXIsIHllYXJfYWRqKV0pCiAgCiAgeWVhcnMgPC0gc2VxKG1pbihzdXJ2ZXlfdW5pdF95ZWFyc1sseWVhcl0pLCBtYXgoc3VydmV5X3VuaXRfeWVhcnNbLHllYXJdKSwgYnkgPSAwLjEpCiAgCiAgeWVhcl9hZGp1c3QgPC0gc2VxKG1pbihzdXJ2ZXlfdW5pdF95ZWFyc1sseWVhcl9hZGpdKSwgbWF4KHN1cnZleV91bml0X3llYXJzWyx5ZWFyX2Fkal0pLCBieSA9IDAuMSkKICAKICB5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdF9hZGRpdGlvbiA8LSBkYXRhLnRhYmxlKHN1cnZleV91bml0ID0gc3VydmV5X3VuaXQubGlzdFtpXSwgeWVhciA9IHllYXJzLCB5ZWFyX2FkaiA9IHllYXJfYWRqdXN0KQogIAogIHllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0IDwtIHJiaW5kKHllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0LCB5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdF9hZGRpdGlvbikKfQoKI2FkZCBjb2xvcnMgYW5kIG5hbWVzIHRvIGZ1bGwgeWVhciBhbmQgc3VydmV5IHVuaXQgY29tYmluYXRpb24gdGFibGUKeWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHQgPC0geWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHRbY29sb3JfbGluaywgb24gPSAic3VydmV5X3VuaXQiXQpgYGAKCgoKR2V0IG1vZGVsIGFzIHByZWRpY3Rpb25zCmBgYHtyfQojZm9yIHBsb3R0aW5nLCBnZXQgbW9kZWwgYXMgcHJlZGljdGlvbnMKYnJheV9jdXJ0aXNfdG90YWxfZ2FtX3ByZWRpY3Rpb25zIDwtIHByZWRpY3QoYnJheV9jdXJ0aXNfdG90YWxfZ2FtLCBzZS5maXQgPSBUUlVFLCBuZXdkYXRhID0geWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHQpCgojbWVyZ2UgaW50byB0YWJsZQp5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdFssYnJheV9nbG1fbW9kX2ZpdCA6PSBicmF5X2N1cnRpc190b3RhbF9nYW1fcHJlZGljdGlvbnMkZml0XVssYnJheV9nbG1fbW9kX2ZpdF9TRSA6PSBicmF5X2N1cnRpc190b3RhbF9nYW1fcHJlZGljdGlvbnMkc2UuZml0XQoKYGBgCgoKUHJvZHVjZSBQbG90IG9mIEdBTSBhbmQgbWVhbiBMTUVSIGxpbmUKYGBge3IgcGxvdCBHQU0gYW5kIG1lYW4gTE1FUiBsaW5lc30KcG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1fbm9sbWVyIDwtIGdncGxvdCgpICsKICMgZ2VvbV9yaWJib24oZGF0YSA9IGxtZXJfYnJheV90b3RhbF9wcmVkaWN0aW9ucywgYWVzKHggPSB5ZWFyLCB5bWluID0gYnJheV9jdXJ0aXNfbG1lcl9wcmVkc19sb3dlckNJLCB5bWF4ID0gYnJheV9jdXJ0aXNfbG1lcl9wcmVkc191cHBlckNJKSwgZmlsbCA9ICJncmV5IiwgYWxwaGEgPSAwLjIpICsKICBnZW9tX3BvaW50KGRhdGEgPSBuYS5vbWl0KGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMucixjb2xzID0gInllYXJfYWRqIiksCiAgICAgICAgICAgICBhZXMoeCA9IHllYXIsCiAgICAgICAgICAgICAgICAgeSA9IGJyYXlfY3VydGlzX2Rpc3NpbWlsYXJpdHlfdG90YWxfbWVhbiwKICAgICAgICAgICAgICAgICBjb2xvciA9IFN1cnZleV9OYW1lX1NlYXNvbiksIGFscGhhID0gMC41LCBzaXplID0gMSkgKwogICAgZ2VvbV9saW5lKGRhdGEgPSBuYS5vbWl0KHllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0LGNvbHMgPSAieWVhcl9hZGoiKSwKICAgICAgICAgICAgIGFlcyh4ID0geWVhciwKICAgICAgICAgICAgICAgICB5ID0gYnJheV9nbG1fbW9kX2ZpdCwKICAgICAgICAgICAgICAgICBjb2xvciA9IFN1cnZleV9OYW1lX1NlYXNvbikpICsKICBnZW9tX3JpYmJvbihkYXRhID0gbmEub21pdCh5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdCxjb2xzID0gInllYXJfYWRqIiksIGFlcyh4ID0geWVhciwgeW1pbj1icmF5X2dsbV9tb2RfZml0LWJyYXlfZ2xtX21vZF9maXRfU0UsIHltYXg9YnJheV9nbG1fbW9kX2ZpdCticmF5X2dsbV9tb2RfZml0X1NFLCBmaWxsID0gIFN1cnZleV9OYW1lX1NlYXNvbiksIGFscGhhPTAuMSkgKyAjYWRkIHN0YW5kYXJkIGVycm9yCiAjIGdlb21fbGluZShkYXRhID0gbG1lcl9icmF5X3RvdGFsX3ByZWRpY3Rpb25zLCBhZXMoeCA9IHllYXIsIHkgPSBicmF5X2N1cnRpc19sbWVyX3ByZWRzKSwgY29sb3IgPSAiYmxhY2siKSArCiAgICBzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzID0gIGNvbG9yX2FscGhhX29yZGVyLCBuYW1lID0gIlN1cnZleSBVbml0IikgKwogIHNjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcyA9ICBjb2xvcl9hbHBoYV9vcmRlciwgZ3VpZGUgPSAibm9uZSIpICsKICB0aGVtZV9jbGFzc2ljKCkgKwogIGxpbXMoeCA9IGMobWluKGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclsseWVhcl0pLG1heChkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLHllYXJdKSksCiAgICAgICB5ID0gYygwLjE1LDAuOSkpICsKICB4bGFiKCJZZWFyIikgKwp5bGFiKCJ0b3RhbCBCQyBkaXNzaW1pbGFyaXR5IikgKwogIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbiA9ICJudWxsIikKCnBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX25vbG1lcgoKZ2dzYXZlKHBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX25vbG1lciwgcGF0aCA9IGhlcmU6OmhlcmUoImZpZ3VyZXMiKSwgZmlsZW5hbWUgPSJwb2ludHNfd2F2eWxpbmVzX2JyYXlfdG90YWxfeWVhcl9yZWR1Y2VkX2dhbV9ub2xtZXIuanBnIiwgaGVpZ2h0ID0gNSwgd2lkdGggPSA1LCB1bml0ID0gImluIikKCiN3aXRoIGxtZXIKCnBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtIDwtIGdncGxvdCgpICsKICBnZW9tX3JpYmJvbihkYXRhID0gbG1lcl9icmF5X3RvdGFsX3ByZWRpY3Rpb25zLCBhZXMoeCA9IHllYXIsIHltaW4gPSBicmF5X2N1cnRpc19sbWVyX3ByZWRzX2xvd2VyQ0ksIHltYXggPSBicmF5X2N1cnRpc19sbWVyX3ByZWRzX3VwcGVyQ0kpLCBmaWxsID0gImdyZXkiLCBhbHBoYSA9IDAuMykgKwogIGdlb21fcG9pbnQoZGF0YSA9IG5hLm9taXQoZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yLCBjb2xzID0gInllYXJfYWRqIiksCiAgICAgICAgICAgICBhZXMoeCA9IHllYXIsCiAgICAgICAgICAgICAgICAgeSA9IGJyYXlfY3VydGlzX2Rpc3NpbWlsYXJpdHlfdG90YWxfbWVhbiwKICAgICAgICAgICAgICAgICBmaWxsID0gU3VydmV5X05hbWVfU2Vhc29uKSwgYWxwaGEgPSAwLjQsIHNpemUgPSAxLCBzaGFwZSA9IDIxLCBjb2xvciA9ICJ3aGl0ZSIpICsKICAgIGdlb21fbGluZShkYXRhID0gbmEub21pdCh5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdCwgY29scyA9ICJ5ZWFyX2FkaiIpLAogICAgICAgICAgICAgYWVzKHggPSB5ZWFyLAogICAgICAgICAgICAgICAgIHkgPSBicmF5X2dsbV9tb2RfZml0LAogICAgICAgICAgICAgICAgIGNvbG9yID0gU3VydmV5X05hbWVfU2Vhc29uKSwgYWxwaGEgPSAwLjYpICsKICBnZW9tX3JpYmJvbihkYXRhID0gbmEub21pdCh5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdCwgY29scyA9ICJ5ZWFyX2FkaiIpLCBhZXMoeCA9IHllYXIsIHltaW49YnJheV9nbG1fbW9kX2ZpdC1icmF5X2dsbV9tb2RfZml0X1NFLCB5bWF4PWJyYXlfZ2xtX21vZF9maXQrYnJheV9nbG1fbW9kX2ZpdF9TRSwgZmlsbCA9ICBTdXJ2ZXlfTmFtZV9TZWFzb24pLCBhbHBoYT0wLjEpICsgI2FkZCBzdGFuZGFyZCBlcnJvcgogIGdlb21fbGluZShkYXRhID0gbG1lcl9icmF5X3RvdGFsX3ByZWRpY3Rpb25zLCBhZXMoeCA9IHllYXIsIHkgPSBicmF5X2N1cnRpc19sbWVyX3ByZWRzKSwgY29sb3IgPSAiYmxhY2siKSArCiAgICBzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzID0gIGNvbG9yX2FscGhhX29yZGVyLCBuYW1lID0gIlN1cnZleSBVbml0IikgKwogIHNjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcyA9ICBjb2xvcl9hbHBoYV9vcmRlciwgZ3VpZGUgPSAibm9uZSIpICsKICB0aGVtZV9jbGFzc2ljKCkgKwogIGxpbXMoeCA9IGMobWluKGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclsseWVhcl0pLG1heChkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLHllYXJdKSkpICsKICB4bGFiKCJZZWFyIikgKwp5bGFiKCLOsi1kaXZlcnNpdHkiKSArCiAgdGhlbWUobGVnZW5kLnBvc2l0aW9uID0gIm51bGwiLCBheGlzLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSwgYXhpcy50aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplID0gMTUpKQoKcG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW0KCmdnc2F2ZShwb2ludHNfd2F2eWxpbmVzX2JyYXlfdG90YWxfeWVhcl9yZWR1Y2VkX2dhbSwgcGF0aCA9IGhlcmU6OmhlcmUoImZpZ3VyZXMiKSwgZmlsZW5hbWUgPSJwb2ludHNfd2F2eWxpbmVzX2JyYXlfdG90YWxfeWVhcl9yZWR1Y2VkX2dhbS5qcGciLCBoZWlnaHQgPSA2LCB3aWR0aCA9IDYsIHVuaXQgPSAiaW4iKQoKI0FMVAojcGxvdCBhbGwsIGJ1dCBzYW1lIGNvbG9yIHNjaGVtZSAoZ3JleSkKcG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1fZ3JleXNjYWxlIDwtIGdncGxvdCgpICsKICBnZW9tX3JpYmJvbihkYXRhID0gbG1lcl9icmF5X3RvdGFsX3ByZWRpY3Rpb25zLCBhZXMoeCA9IHllYXIsIHltaW4gPSBicmF5X2N1cnRpc19sbWVyX3ByZWRzX2xvd2VyQ0ksIHltYXggPSBicmF5X2N1cnRpc19sbWVyX3ByZWRzX3VwcGVyQ0kpLCBmaWxsID0gImdyZXkiLCBhbHBoYSA9IDAuMykgKwogIGdlb21fcG9pbnQoZGF0YSA9IGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuciwKICAgICAgICAgICAgIGFlcyh4ID0geWVhciwKICAgICAgICAgICAgICAgICB5ID0gYnJheV9jdXJ0aXNfZGlzc2ltaWxhcml0eV90b3RhbF9tZWFuLAogICAgICAgICAgICAgICAgIGZpbGwgPSBTdXJ2ZXlfTmFtZV9TZWFzb24pLCBhbHBoYSA9IDAuNCwgc2l6ZSA9IDEsIHNoYXBlID0gMjEsIGNvbG9yID0gIndoaXRlIikgKwogICAgZ2VvbV9saW5lKGRhdGEgPSB5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdCwKICAgICAgICAgICAgIGFlcyh4ID0geWVhciwKICAgICAgICAgICAgICAgICB5ID0gYnJheV9nbG1fbW9kX2ZpdCwKICAgICAgICAgICAgICAgICBjb2xvciA9IFN1cnZleV9OYW1lX1NlYXNvbiksIGFscGhhID0gMC42KSArCiAgZ2VvbV9yaWJib24oZGF0YSA9IHllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0LCBhZXMoeCA9IHllYXIsIHltaW49YnJheV9nbG1fbW9kX2ZpdC1icmF5X2dsbV9tb2RfZml0X1NFLCB5bWF4PWJyYXlfZ2xtX21vZF9maXQrYnJheV9nbG1fbW9kX2ZpdF9TRSwgZmlsbCA9ICBTdXJ2ZXlfTmFtZV9TZWFzb24pLCBhbHBoYT0wLjEpICsgI2FkZCBzdGFuZGFyZCBlcnJvcgogIGdlb21fbGluZShkYXRhID0gbG1lcl9icmF5X3RvdGFsX3ByZWRpY3Rpb25zLCBhZXMoeCA9IHllYXIsIHkgPSBicmF5X2N1cnRpc19sbWVyX3ByZWRzKSwgY29sb3IgPSAiYmxhY2siKSArCiAgICBzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzID0gIHJlcCgiYmxhY2siLCB0aW1lcyA9IGxlbmd0aCh1bmlxdWUoZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yJFN1cnZleV9OYW1lX1NlYXNvbikpKSwgbmFtZSA9ICJTdXJ2ZXkgVW5pdCIpICsKICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSAgcmVwKCJibGFjayIsIHRpbWVzID0gbGVuZ3RoKHVuaXF1ZShkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnIkU3VydmV5X05hbWVfU2Vhc29uKSkpLCBndWlkZSA9ICJub25lIikgKwogIHRoZW1lX2NsYXNzaWMoKSArCiAgbGltcyh4ID0gYyhtaW4oZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyx5ZWFyXSksbWF4KGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclsseWVhcl0pKQogICAgICAgKSArCiAgeGxhYigiWWVhciIpICsKeWxhYigizrItZGl2ZXJzaXR5IikgKwogIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbiA9ICJudWxsIiwgYXhpcy50ZXh0ID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNSksIGF4aXMudGl0bGUgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSkKCnBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX2dyZXlzY2FsZQoKZ2dzYXZlKHBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX2dyZXlzY2FsZSwgcGF0aCA9IGhlcmU6OmhlcmUoImZpZ3VyZXMiKSwgZmlsZW5hbWUgPSJwb2ludHNfd2F2eWxpbmVzX2JyYXlfdG90YWxfeWVhcl9yZWR1Y2VkX2dhbV9ncmV5c2NhbGUuanBnIiwgaGVpZ2h0ID0gNiwgd2lkdGggPSA2LCB1bml0ID0gImluIikKYGBgCkFsdGVybmF0aXZlLCBjb2xvciBieSB0cmVuZApgYGB7ciBjb2xvciB3YXZ5IGxpbmVzIGJ5IHRyZW5kfQoKcG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1fY29sb3JieXRyZW5kIDwtIGdncGxvdCgpICsKICBnZW9tX3JpYmJvbihkYXRhID0gbG1lcl9icmF5X3RvdGFsX3ByZWRpY3Rpb25zLCBhZXMoeCA9IHllYXIsIHltaW4gPSBicmF5X2N1cnRpc19sbWVyX3ByZWRzX2xvd2VyQ0ksIHltYXggPSBicmF5X2N1cnRpc19sbWVyX3ByZWRzX3VwcGVyQ0kpLCBmaWxsID0gImdyZXkiLCBhbHBoYSA9IDAuMykgKwogIGdlb21fcG9pbnQoZGF0YSA9IG5hLm9taXQoZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yLCBjb2xzID0gInllYXJfYWRqIiksCiAgICAgICAgICAgICBhZXMoeCA9IHllYXIsCiAgICAgICAgICAgICAgICAgeSA9IGJyYXlfY3VydGlzX2Rpc3NpbWlsYXJpdHlfdG90YWxfbWVhbiwKICAgICAgICAgICAgICAgICBmaWxsID0gU3VydmV5X05hbWVfU2Vhc29uKSwgYWxwaGEgPSAwLjQsIHNpemUgPSAxLCBzaGFwZSA9IDIxLCBjb2xvciA9ICJ3aGl0ZSIpICsKICAgIGdlb21fbGluZShkYXRhID0gbmEub21pdCh5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdCwgY29scyA9ICJ5ZWFyX2FkaiIpLAogICAgICAgICAgICAgYWVzKHggPSB5ZWFyLAogICAgICAgICAgICAgICAgIHkgPSBicmF5X2dsbV9tb2RfZml0LAogICAgICAgICAgICAgICAgIGNvbG9yID0gU3VydmV5X05hbWVfU2Vhc29uKSwgYWxwaGEgPSAwLjYpICsKICBnZW9tX3JpYmJvbihkYXRhID0gbmEub21pdCh5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdCwgY29scyA9ICJ5ZWFyX2FkaiIpLCBhZXMoeCA9IHllYXIsIHltaW49YnJheV9nbG1fbW9kX2ZpdC1icmF5X2dsbV9tb2RfZml0X1NFLCB5bWF4PWJyYXlfZ2xtX21vZF9maXQrYnJheV9nbG1fbW9kX2ZpdF9TRSwgZmlsbCA9ICBTdXJ2ZXlfTmFtZV9TZWFzb24pLCBhbHBoYT0wLjEpICsgI2FkZCBzdGFuZGFyZCBlcnJvcgogIGdlb21fbGluZShkYXRhID0gbG1lcl9icmF5X3RvdGFsX3ByZWRpY3Rpb25zLCBhZXMoeCA9IHllYXIsIHkgPSBicmF5X2N1cnRpc19sbWVyX3ByZWRzKSwgY29sb3IgPSAiYmxhY2siKSArCiAgICBzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzID0gIGNvbG9yX2FscGhhX29yZGVyX2J5dHJlbmQsIG5hbWUgPSAiU3VydmV5IFVuaXQiKSArCiAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gIGNvbG9yX2FscGhhX29yZGVyX2J5dHJlbmQsIGd1aWRlID0gIm5vbmUiKSArCiAgdGhlbWVfY2xhc3NpYygpICsKICBsaW1zKHggPSBjKG1pbihkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLHllYXJdKSxtYXgoZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyx5ZWFyXSkpLHkgPSBjKDAuNDksMC45NSkpICsKICB4bGFiKCJZZWFyIikgKwp5bGFiKCLOsi1kaXZlcnNpdHkiKSArCiAgdGhlbWUobGVnZW5kLnBvc2l0aW9uID0gIm51bGwiLCBheGlzLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSwgYXhpcy50aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplID0gMTUpKQoKcG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1fY29sb3JieXRyZW5kCgpnZ3NhdmUocG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1fY29sb3JieXRyZW5kLCBwYXRoID0gaGVyZTo6aGVyZSgiZmlndXJlcyIpLCBmaWxlbmFtZSA9InBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX2NvbG9yYnl0cmVuZC5qcGciLCBoZWlnaHQgPSA2LCB3aWR0aCA9IDExLCB1bml0ID0gImluIikKCgpwb2ludHNfd2F2eWxpbmVzX2JyYXlfdG90YWxfeWVhcl9yZWR1Y2VkX2dhbV9jb2xvcmJ5dHJlbmRfc3F1aWdnbGVzb25seSA8LSBnZ3Bsb3QoKSArCiAgZ2VvbV9wb2ludChkYXRhID0gbmEub21pdChkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnIsIGNvbHMgPSAieWVhcl9hZGoiKSwKICAgICAgICAgICAgIGFlcyh4ID0geWVhciwKICAgICAgICAgICAgICAgICB5ID0gYnJheV9jdXJ0aXNfZGlzc2ltaWxhcml0eV90b3RhbF9tZWFuLAogICAgICAgICAgICAgICAgIGZpbGwgPSBTdXJ2ZXlfTmFtZV9TZWFzb24pLCBhbHBoYSA9IDAuNCwgc2l6ZSA9IDEsIHNoYXBlID0gMjEsIGNvbG9yID0gIndoaXRlIikgKwogICAgZ2VvbV9saW5lKGRhdGEgPSBuYS5vbWl0KHllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0LCBjb2xzID0gInllYXJfYWRqIiksCiAgICAgICAgICAgICBhZXMoeCA9IHllYXIsCiAgICAgICAgICAgICAgICAgeSA9IGJyYXlfZ2xtX21vZF9maXQsCiAgICAgICAgICAgICAgICAgY29sb3IgPSBTdXJ2ZXlfTmFtZV9TZWFzb24pLCBhbHBoYSA9IDAuNikgKwogIGdlb21fcmliYm9uKGRhdGEgPSBuYS5vbWl0KHllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0LCBjb2xzID0gInllYXJfYWRqIiksIGFlcyh4ID0geWVhciwgeW1pbj1icmF5X2dsbV9tb2RfZml0LWJyYXlfZ2xtX21vZF9maXRfU0UsIHltYXg9YnJheV9nbG1fbW9kX2ZpdCticmF5X2dsbV9tb2RfZml0X1NFLCBmaWxsID0gIFN1cnZleV9OYW1lX1NlYXNvbiksIGFscGhhPTAuMSkgKyAjYWRkIHN0YW5kYXJkIGVycm9yCiAgICBzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzID0gIGNvbG9yX2FscGhhX29yZGVyX2J5dHJlbmQsIG5hbWUgPSAiU3VydmV5IFVuaXQiKSArCiAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gIGNvbG9yX2FscGhhX29yZGVyX2J5dHJlbmQsIGd1aWRlID0gIm5vbmUiKSArCiAgdGhlbWVfY2xhc3NpYygpICsKICBsaW1zKHggPSBjKG1pbihkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLHllYXJdKSxtYXgoZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyx5ZWFyXSkpLCB5ID0gYygwLjQ5LDAuOTUpKSArCiAgeGxhYigiWWVhciIpICsKeWxhYigizrItZGl2ZXJzaXR5IikgKwogIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbiA9ICJudWxsIiwgYXhpcy50ZXh0ID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNSksIGF4aXMudGl0bGUgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSkKCnBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX2NvbG9yYnl0cmVuZF9zcXVpZ2dsZXNvbmx5CgpnZ3NhdmUocG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1fY29sb3JieXRyZW5kX3NxdWlnZ2xlc29ubHksIHBhdGggPSBoZXJlOjpoZXJlKCJmaWd1cmVzIiksIGZpbGVuYW1lID0icG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1fY29sb3JieXRyZW5kX3NxdWlnZ2xlc29ubHkuanBnIiwgaGVpZ2h0ID0gNiwgd2lkdGggPSAxMSwgdW5pdCA9ICJpbiIpCmBgYAoKYGBge3IgcGxvdCBlYWNoIGluZGVwZW5kZW50bHl9CiNwbG90IGVhY2ggaW5kZXBlbmRlbnRseSBmb3Igc3VwcGxlbWVudAojYWxsIHN1cnZleSBuYW1lcyA9IAphbGxfc3VydmV5X25hbWVzIDwtIHNvcnQodW5pcXVlKGNvbG9yX2xpbmskU3VydmV5X05hbWVfU2Vhc29uKSkKI2xpc3Qgb2YgcGxvdHMKcG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1faW5kaXZpZHVhbCA8LSBsaXN0KCkKZm9yIChpIGluIDE6bGVuZ3RoKGFsbF9zdXJ2ZXlfbmFtZXMpKSB7CnBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX2luZGl2aWR1YWxbW2ldXSA8LSBnZ3Bsb3QoKSArCiAgZ2VvbV9wb2ludChkYXRhID0gZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yW1N1cnZleV9OYW1lX1NlYXNvbiA9PSBhbGxfc3VydmV5X25hbWVzW2ldXSwKICAgICAgICAgICAgIGFlcyh4ID0geWVhciwKICAgICAgICAgICAgICAgICB5ID0gYnJheV9jdXJ0aXNfZGlzc2ltaWxhcml0eV90b3RhbF9tZWFuKSwgYWxwaGEgPSAwLjQsIGNvbG9yID0gImJsYWNrIikgKwogICAgZ2VvbV9saW5lKGRhdGEgPSB5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdFtTdXJ2ZXlfTmFtZV9TZWFzb24gPT0gYWxsX3N1cnZleV9uYW1lc1tpXV0sCiAgICAgICAgICAgICBhZXMoeCA9IHllYXIsCiAgICAgICAgICAgICAgICAgeSA9IGJyYXlfZ2xtX21vZF9maXQpLCBhbHBoYSA9IDAuNikgKwogIGdlb21fcmliYm9uKGRhdGEgPSB5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdFtTdXJ2ZXlfTmFtZV9TZWFzb24gPT0gYWxsX3N1cnZleV9uYW1lc1tpXV0sIGFlcyh4ID0geWVhciwgeW1pbj1icmF5X2dsbV9tb2RfZml0LWJyYXlfZ2xtX21vZF9maXRfU0UsIHltYXg9YnJheV9nbG1fbW9kX2ZpdCticmF5X2dsbV9tb2RfZml0X1NFKSwgYWxwaGE9MC4xKSArICNhZGQgc3RhbmRhcmQgZXJyb3IKICB0aGVtZV9jbGFzc2ljKCkgKwojICBsaW1zKHggPSBjKG1pbihkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbU3VydmV5X05hbWVfU2Vhc29uID09IGFsbF9zdXJ2ZXlfbmFtZXNbaV0seWVhcl0pLG1heChkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbU3VydmV5X05hbWVfU2Vhc29uID09IGFsbF9zdXJ2ZXlfbmFtZXNbaV0seWVhcl0pKSwKIyAgICAgICB5ID0gYygwLjE1LDAuOSkpICsKICB4bGFiKCJZZWFyIikgKwp5bGFiKCJiZXRhLWRpdmVyc2l0eSIpICsKICBmYWNldF93cmFwKH5TdXJ2ZXlfTmFtZV9TZWFzb24sIG5jb2wgPSA1KSArCiAgdGhlbWUobGVnZW5kLnBvc2l0aW9uID0gIm51bGwiLCBheGlzLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSwgYXhpcy50aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplID0gMTUpKQoKcHJpbnQocG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1faW5kaXZpZHVhbFtbaV1dKQoKfQpzYXZlUkRTKHBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX2luZGl2aWR1YWwsIGhlcmU6OmhlcmUoImZpZ3VyZXMiLCJwb2ludHNfd2F2eWxpbmVzX2JyYXlfdG90YWxfeWVhcl9yZWR1Y2VkX2dhbV9pbmRpdmlkdWFsLlJkcyIpKQoKI3ByaW50IHRvIHBkZgpsaWJyYXJ5KGdyaWRFeHRyYSkKCmdnc2F2ZSgKICAgZmlsZW5hbWUgPSBoZXJlOjpoZXJlKCJmaWd1cmVzIiwicG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1faW5kaXZpZHVhbC5wZGYiKSwgCiAgIHBsb3QgPSBtYXJyYW5nZUdyb2IocG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1faW5kaXZpZHVhbCwgbnJvdz0xLCBuY29sPTEpLCAKICAgd2lkdGggPSA4LjUsIGhlaWdodCA9IDExCikKCiNBbHRlcm5hdGl2ZWx5LCBzcGxpdCBpbnRvIDIgYW5kIHVzZSBmYWNldAojZmlyc3QgMjQKcG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1faW5kaXZpZHVhbF9mYWNldF8xXzI0IDwtIGdncGxvdCgpICsKICBnZW9tX3BvaW50KGRhdGEgPSBkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbU3VydmV5X05hbWVfU2Vhc29uICAlaW4lIGFsbF9zdXJ2ZXlfbmFtZXNbYygxOjI0KV1dLAogICAgICAgICAgICAgYWVzKHggPSB5ZWFyLAogICAgICAgICAgICAgICAgIHkgPSBicmF5X2N1cnRpc19kaXNzaW1pbGFyaXR5X3RvdGFsX21lYW4pLCBhbHBoYSA9IDAuNykgKwogICAgZ2VvbV9saW5lKGRhdGEgPSB5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdFtTdXJ2ZXlfTmFtZV9TZWFzb24gICAlaW4lIGFsbF9zdXJ2ZXlfbmFtZXNbYygxOjI0KV1dLAogICAgICAgICAgICAgYWVzKHggPSB5ZWFyLAogICAgICAgICAgICAgICAgIHkgPSBicmF5X2dsbV9tb2RfZml0KSkgKwogIGdlb21fcmliYm9uKGRhdGEgPSB5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdFtTdXJ2ZXlfTmFtZV9TZWFzb24gICVpbiUgYWxsX3N1cnZleV9uYW1lc1tjKDE6MjQpXV0sCmFlcyh4ID0geWVhciwgeW1pbj1icmF5X2dsbV9tb2RfZml0LWJyYXlfZ2xtX21vZF9maXRfU0UsIHltYXg9YnJheV9nbG1fbW9kX2ZpdCticmF5X2dsbV9tb2RfZml0X1NFKSwgYWxwaGE9MC41KSArICNhZGQgc3RhbmRhcmQgZXJyb3IKICB0aGVtZV9jbGFzc2ljKCkgKwogIHhsYWIoIlllYXIiKSArCnlsYWIoIs6yLWRpdmVyc2l0eSIpICsKICBmYWNldF93cmFwKH5TdXJ2ZXlfTmFtZV9TZWFzb24sIG5jb2wgPSA0LCBzY2FsZXMgPSAiZnJlZSIpICsKICB0aGVtZShheGlzLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDgpLCBheGlzLnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemUgPSAxMikpCgpnZ3NhdmUocG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1faW5kaXZpZHVhbF9mYWNldF8xXzI0LCBwYXRoID0gIGhlcmU6OmhlcmUoImZpZ3VyZXMiKSwgZmlsZW5hbWUgPSAicG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1faW5kaXZpZHVhbF9mYWNldF8xXzI0LnBuZyIsIGhlaWdodCA9IDExLCB3aWR0aCA9IDkpCgpwb2ludHNfd2F2eWxpbmVzX2JyYXlfdG90YWxfeWVhcl9yZWR1Y2VkX2dhbV9pbmRpdmlkdWFsX2ZhY2V0XzI1XzQyIDwtIGdncGxvdCgpICsKICBnZW9tX3BvaW50KGRhdGEgPSBkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbU3VydmV5X05hbWVfU2Vhc29uICAlaW4lIGFsbF9zdXJ2ZXlfbmFtZXNbYygyNTo0MildXSwKICAgICAgICAgICAgIGFlcyh4ID0geWVhciwKICAgICAgICAgICAgICAgICB5ID0gYnJheV9jdXJ0aXNfZGlzc2ltaWxhcml0eV90b3RhbF9tZWFuKSwgYWxwaGEgPSAwLjcpICsKICAgIGdlb21fbGluZShkYXRhID0geWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHRbU3VydmV5X05hbWVfU2Vhc29uICAgJWluJSBhbGxfc3VydmV5X25hbWVzW2MoMjU6NDIpXV0sCiAgICAgICAgICAgICBhZXMoeCA9IHllYXIsCiAgICAgICAgICAgICAgICAgeSA9IGJyYXlfZ2xtX21vZF9maXQpKSArCiAgZ2VvbV9yaWJib24oZGF0YSA9IHllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0W1N1cnZleV9OYW1lX1NlYXNvbiAgJWluJSBhbGxfc3VydmV5X25hbWVzW2MoMjU6NDIpXV0sCmFlcyh4ID0geWVhciwgeW1pbj1icmF5X2dsbV9tb2RfZml0LWJyYXlfZ2xtX21vZF9maXRfU0UsIHltYXg9YnJheV9nbG1fbW9kX2ZpdCticmF5X2dsbV9tb2RfZml0X1NFKSwgYWxwaGE9MC41KSArICNhZGQgc3RhbmRhcmQgZXJyb3IKICB0aGVtZV9jbGFzc2ljKCkgKwogIHhsYWIoIlllYXIiKSArCnlsYWIoIs6yLWRpdmVyc2l0eSIpICsKICBmYWNldF93cmFwKH5TdXJ2ZXlfTmFtZV9TZWFzb24sIG5jb2wgPSA0LCBzY2FsZXMgPSAiZnJlZSIpICsKICB0aGVtZShheGlzLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDgpLCBheGlzLnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemUgPSAxMikpCgpnZ3NhdmUocG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1faW5kaXZpZHVhbF9mYWNldF8yNV80MiwgcGF0aCA9ICBoZXJlOjpoZXJlKCJmaWd1cmVzIiksIGZpbGVuYW1lID0gInBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX2luZGl2aWR1YWxfZmFjZXRfMjVfNDIucG5nIiwgaGVpZ2h0ID0gMTEsIHdpZHRoID0gOSkKCgpgYGAKClJlZ2lvbiBieSByZWdpb24gYnVpbGQgdXAKYGBge3IgcmVnaW9uIGJ1aWxkIHVwfQojRUJTIE9OTFkgKHNlZ3VlIGZvciBjaGFwdGVyIDMpCnBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX2NvbG9yYnl0cmVuZF9FQlNvbmx5IDwtIGdncGxvdCgpICsKICBnZW9tX3BvaW50KGRhdGEgPSBuYS5vbWl0KGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMucltzdXJ2ZXlfdW5pdCA9PSAiRUJTIl0sIGNvbHMgPSAieWVhcl9hZGoiKSwKICAgICAgICAgICAgIGFlcyh4ID0geWVhciwKICAgICAgICAgICAgICAgICB5ID0gYnJheV9jdXJ0aXNfZGlzc2ltaWxhcml0eV90b3RhbF9tZWFuLAogICAgICAgICAgICAgICAgIGZpbGwgPSBTdXJ2ZXlfTmFtZV9TZWFzb24pLCBzaXplID0gMiwgc2hhcGUgPSAyMSwgY29sb3IgPSAid2hpdGUiKSArCiAgICBnZW9tX2xpbmUoZGF0YSA9IG5hLm9taXQoeWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHRbc3VydmV5X3VuaXQgPT0gIkVCUyJdLCBjb2xzID0gInllYXJfYWRqIiksCiAgICAgICAgICAgICBhZXMoeCA9IHllYXIsCiAgICAgICAgICAgICAgICAgeSA9IGJyYXlfZ2xtX21vZF9maXQsCiAgICAgICAgICAgICAgICAgY29sb3IgPSBTdXJ2ZXlfTmFtZV9TZWFzb24pKSArCiAgZ2VvbV9yaWJib24oZGF0YSA9IG5hLm9taXQoeWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHRbc3VydmV5X3VuaXQgPT0gIkVCUyJdLCBjb2xzID0gInllYXJfYWRqIiksIGFlcyh4ID0geWVhciwgeW1pbj1icmF5X2dsbV9tb2RfZml0LWJyYXlfZ2xtX21vZF9maXRfU0UsIHltYXg9YnJheV9nbG1fbW9kX2ZpdCticmF5X2dsbV9tb2RfZml0X1NFLCBmaWxsID0gIFN1cnZleV9OYW1lX1NlYXNvbiksIGFscGhhPTAuMSkgKyAjYWRkIHN0YW5kYXJkIGVycm9yCiAgICBzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzID0gICIjY2JiZmRlIiwgbmFtZSA9ICJTdXJ2ZXkgVW5pdCIpICsKICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSAgIiNjYmJmZGUiLCBndWlkZSA9ICJub25lIikgKwogIHRoZW1lX2NsYXNzaWMoKSArCiAgbGltcyh4ID0gYyhtaW4oZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyx5ZWFyXSksbWF4KGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclsseWVhcl0pKSwgeSA9IGMoMC40OSwwLjk1KSkgKwogIHhsYWIoIlllYXIiKSArCnlsYWIoIs6yLWRpdmVyc2l0eSIpICsKICB0aGVtZShsZWdlbmQucG9zaXRpb24gPSAibnVsbCIsIGF4aXMudGV4dCA9IGVsZW1lbnRfdGV4dChzaXplID0gMTUpLCBheGlzLnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNSkpCgpwb2ludHNfd2F2eWxpbmVzX2JyYXlfdG90YWxfeWVhcl9yZWR1Y2VkX2dhbV9jb2xvcmJ5dHJlbmRfRUJTb25seQoKZ2dzYXZlKHBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX2NvbG9yYnl0cmVuZF9FQlNvbmx5LCBwYXRoID0gaGVyZTo6aGVyZSgiZmlndXJlcyIpLCBmaWxlbmFtZSA9InBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX2NvbG9yYnl0cmVuZF9FQlNvbmx5LmpwZyIsIGhlaWdodCA9IDYsIHdpZHRoID0gMTEsIHVuaXQgPSAiaW4iKQoKCiNqdXN0IG9uZSByZWdpb24gdGhhdCdzIGhvbW9nZW5pemluZyAoU0VVUy1zdW1tZXIpCnBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX2NvbG9yYnl0cmVuZF9zZXVzX3N1bW1lciA8LSBnZ3Bsb3QoKSArCiAgZ2VvbV9wb2ludChkYXRhID0gbmEub21pdChkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbc3VydmV5X3VuaXQgPT0gIlNFVVMtc3VtbWVyIl0sIGNvbHMgPSAieWVhcl9hZGoiKSwKICAgICAgICAgICAgIGFlcyh4ID0geWVhciwKICAgICAgICAgICAgICAgICB5ID0gYnJheV9jdXJ0aXNfZGlzc2ltaWxhcml0eV90b3RhbF9tZWFuLAogICAgICAgICAgICAgICAgIGZpbGwgPSBTdXJ2ZXlfTmFtZV9TZWFzb24pLCBzaXplID0gMiwgc2hhcGUgPSAyMSwgY29sb3IgPSAid2hpdGUiKSArCiAgICBnZW9tX2xpbmUoZGF0YSA9IG5hLm9taXQoeWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHRbc3VydmV5X3VuaXQgPT0gIlNFVVMtc3VtbWVyIl0sIGNvbHMgPSAieWVhcl9hZGoiKSwKICAgICAgICAgICAgIGFlcyh4ID0geWVhciwKICAgICAgICAgICAgICAgICB5ID0gYnJheV9nbG1fbW9kX2ZpdCwKICAgICAgICAgICAgICAgICBjb2xvciA9IFN1cnZleV9OYW1lX1NlYXNvbikpICsKICBnZW9tX3JpYmJvbihkYXRhID0gbmEub21pdCh5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdFtzdXJ2ZXlfdW5pdCA9PSAiU0VVUy1zdW1tZXIiXSwgY29scyA9ICJ5ZWFyX2FkaiIpLCBhZXMoeCA9IHllYXIsIHltaW49YnJheV9nbG1fbW9kX2ZpdC1icmF5X2dsbV9tb2RfZml0X1NFLCB5bWF4PWJyYXlfZ2xtX21vZF9maXQrYnJheV9nbG1fbW9kX2ZpdF9TRSwgZmlsbCA9ICBTdXJ2ZXlfTmFtZV9TZWFzb24pLCBhbHBoYT0wLjEpICsgI2FkZCBzdGFuZGFyZCBlcnJvcgogICAgc2NhbGVfY29sb3JfbWFudWFsKHZhbHVlcyA9ICAiI2U3YWM1YiIsIG5hbWUgPSAiU3VydmV5IFVuaXQiKSArCiAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gICIjZTdhYzViIiwgZ3VpZGUgPSAibm9uZSIpICsKICB0aGVtZV9jbGFzc2ljKCkgKwogIGxpbXMoeCA9IGMobWluKGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclsseWVhcl0pLG1heChkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLHllYXJdKSksIHkgPSBjKDAuNDksMC45NSkpICsKICB4bGFiKCJZZWFyIikgKwp5bGFiKCLOsi1kaXZlcnNpdHkiKSArCiAgdGhlbWUobGVnZW5kLnBvc2l0aW9uID0gIm51bGwiLCBheGlzLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSwgYXhpcy50aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplID0gMTUpKQoKcG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1fY29sb3JieXRyZW5kX3NldXNfc3VtbWVyCgpnZ3NhdmUocG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1fY29sb3JieXRyZW5kX3NldXNfc3VtbWVyLCBwYXRoID0gaGVyZTo6aGVyZSgiZmlndXJlcyIpLCBmaWxlbmFtZSA9InBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX2NvbG9yYnl0cmVuZF9zZXVzX3N1bW1lci5qcGciLCBoZWlnaHQgPSA2LCB3aWR0aCA9IDExLCB1bml0ID0gImluIikKCgoKI2p1c3Qgb25lIHJlZ2lvbiB0aGF0J3MgZGlmZmVyZW50aWF0aW5nIGFuZCBvbmUgdGhhdCdzIGhvbW9nZW5pemluZwpwb2ludHNfd2F2eWxpbmVzX2JyYXlfdG90YWxfeWVhcl9yZWR1Y2VkX2dhbV9jb2xvcmJ5dHJlbmRfTlpfV0NTSV9TRVVTX3N1bW1lcm9ubHkgPC0gZ2dwbG90KCkgKwogIGdlb21fcG9pbnQoZGF0YSA9IG5hLm9taXQoZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yW3N1cnZleV91bml0ICVpbiUgYygiTlotV0NTSSIsIlNFVVMtc3VtbWVyIildLCBjb2xzID0gInllYXJfYWRqIiksCiAgICAgICAgICAgICBhZXMoeCA9IHllYXIsCiAgICAgICAgICAgICAgICAgeSA9IGJyYXlfY3VydGlzX2Rpc3NpbWlsYXJpdHlfdG90YWxfbWVhbiwKICAgICAgICAgICAgICAgICBmaWxsID0gU3VydmV5X05hbWVfU2Vhc29uKSwgc2l6ZSA9IDIsIHNoYXBlID0gMjEsIGNvbG9yID0gIndoaXRlIikgKwogICAgZ2VvbV9saW5lKGRhdGEgPSBuYS5vbWl0KHllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0W3N1cnZleV91bml0ICVpbiUgYygiTlotV0NTSSIsIlNFVVMtc3VtbWVyIildLCBjb2xzID0gInllYXJfYWRqIiksCiAgICAgICAgICAgICBhZXMoeCA9IHllYXIsCiAgICAgICAgICAgICAgICAgeSA9IGJyYXlfZ2xtX21vZF9maXQsCiAgICAgICAgICAgICAgICAgY29sb3IgPSBTdXJ2ZXlfTmFtZV9TZWFzb24pKSArCiAgZ2VvbV9yaWJib24oZGF0YSA9IG5hLm9taXQoeWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHRbc3VydmV5X3VuaXQgJWluJSBjKCJOWi1XQ1NJIiwiU0VVUy1zdW1tZXIiKV0sIGNvbHMgPSAieWVhcl9hZGoiKSwgYWVzKHggPSB5ZWFyLCB5bWluPWJyYXlfZ2xtX21vZF9maXQtYnJheV9nbG1fbW9kX2ZpdF9TRSwgeW1heD1icmF5X2dsbV9tb2RfZml0K2JyYXlfZ2xtX21vZF9maXRfU0UsIGZpbGwgPSAgU3VydmV5X05hbWVfU2Vhc29uKSwgYWxwaGE9MC4xKSArICNhZGQgc3RhbmRhcmQgZXJyb3IKICAgIHNjYWxlX2NvbG9yX21hbnVhbCh2YWx1ZXMgPSAgYygiI2U3YWM1YiIsIiM5MWM4NzQiKSwgbmFtZSA9ICJTdXJ2ZXkgVW5pdCIpICsKICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSAgYygiI2U3YWM1YiIsIiM5MWM4NzQiKSwgZ3VpZGUgPSAibm9uZSIpICsKICB0aGVtZV9jbGFzc2ljKCkgKwogIGxpbXMoeCA9IGMobWluKGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclsseWVhcl0pLG1heChkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLHllYXJdKSksIHkgPSBjKDAuNDksMC45NSkpICsKICB4bGFiKCJZZWFyIikgKwp5bGFiKCLOsi1kaXZlcnNpdHkiKSArCiAgdGhlbWUobGVnZW5kLnBvc2l0aW9uID0gIm51bGwiLCBheGlzLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSwgYXhpcy50aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplID0gMTUpKQoKcG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1fY29sb3JieXRyZW5kX05aX1dDU0lfU0VVU19zdW1tZXJvbmx5CgpnZ3NhdmUocG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1fY29sb3JieXRyZW5kX05aX1dDU0lfU0VVU19zdW1tZXJvbmx5LCBwYXRoID0gaGVyZTo6aGVyZSgiZmlndXJlcyIpLCBmaWxlbmFtZSA9InBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX2NvbG9yYnl0cmVuZF9OWl9XQ1NJX1NFVVNfc3VtbWVyb25seS5qcGciLCBoZWlnaHQgPSA2LCB3aWR0aCA9IDExLCB1bml0ID0gImluIikKCiNqdXN0IHRocmVlIHJlZ2lvbnMsIGFub3RoZXIgb25lIHRoYXQncyBzdGFibGUsIChORVVTLUZhbGwpCnBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX2NvbG9yYnl0cmVuZF9OWl9XQ1NJX1NFVVNfc3VtbWVyX05FVVNmYWxsb25seSA8LSBnZ3Bsb3QoKSArCiAgZ2VvbV9wb2ludChkYXRhID0gbmEub21pdChkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbc3VydmV5X3VuaXQgJWluJSBjKCJOWi1XQ1NJIiwiU0VVUy1zdW1tZXIiLCAiTkVVUy1GYWxsIildLCBjb2xzID0gInllYXJfYWRqIiksCiAgICAgICAgICAgICBhZXMoeCA9IHllYXIsCiAgICAgICAgICAgICAgICAgeSA9IGJyYXlfY3VydGlzX2Rpc3NpbWlsYXJpdHlfdG90YWxfbWVhbiwKICAgICAgICAgICAgICAgICBmaWxsID0gU3VydmV5X05hbWVfU2Vhc29uKSwgc2l6ZSA9IDIsIHNoYXBlID0gMjEsIGNvbG9yID0gIndoaXRlIikgKwogICAgZ2VvbV9saW5lKGRhdGEgPSBuYS5vbWl0KHllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0W3N1cnZleV91bml0ICVpbiUgYygiTlotV0NTSSIsIlNFVVMtc3VtbWVyIiwgIk5FVVMtRmFsbCIpXSwgY29scyA9ICJ5ZWFyX2FkaiIpLAogICAgICAgICAgICAgYWVzKHggPSB5ZWFyLAogICAgICAgICAgICAgICAgIHkgPSBicmF5X2dsbV9tb2RfZml0LAogICAgICAgICAgICAgICAgIGNvbG9yID0gU3VydmV5X05hbWVfU2Vhc29uKSkgKwogIGdlb21fcmliYm9uKGRhdGEgPSBuYS5vbWl0KHllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0W3N1cnZleV91bml0ICVpbiUgYygiTlotV0NTSSIsIlNFVVMtc3VtbWVyIiwgIk5FVVMtRmFsbCIpXSwgY29scyA9ICJ5ZWFyX2FkaiIpLCBhZXMoeCA9IHllYXIsIHltaW49YnJheV9nbG1fbW9kX2ZpdC1icmF5X2dsbV9tb2RfZml0X1NFLCB5bWF4PWJyYXlfZ2xtX21vZF9maXQrYnJheV9nbG1fbW9kX2ZpdF9TRSwgZmlsbCA9ICBTdXJ2ZXlfTmFtZV9TZWFzb24pLCBhbHBoYT0wLjEpICsgI2FkZCBzdGFuZGFyZCBlcnJvcgogICAgc2NhbGVfY29sb3JfbWFudWFsKHZhbHVlcyA9ICBjKCIjY2JiZmRlIiwiI2U3YWM1YiIsICIjOTFjODc0IiksIG5hbWUgPSAiU3VydmV5IFVuaXQiKSArCiAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gIGMoIiNjYmJmZGUiLCIjZTdhYzViIiwgIiM5MWM4NzQiKSwgZ3VpZGUgPSAibm9uZSIpICsKICB0aGVtZV9jbGFzc2ljKCkgKwogIGxpbXMoeCA9IGMobWluKGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclsseWVhcl0pLG1heChkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLHllYXJdKSksIHkgPSBjKDAuNDksMC45NSkpICsKICB4bGFiKCJZZWFyIikgKwp5bGFiKCLOsi1kaXZlcnNpdHkiKSArCiAgdGhlbWUobGVnZW5kLnBvc2l0aW9uID0gIm51bGwiLCBheGlzLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSwgYXhpcy50aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplID0gMTUpKQoKcG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1fY29sb3JieXRyZW5kX05aX1dDU0lfU0VVU19zdW1tZXJfTkVVU2ZhbGxvbmx5CgpnZ3NhdmUocG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1fY29sb3JieXRyZW5kX05aX1dDU0lfU0VVU19zdW1tZXJfTkVVU2ZhbGxvbmx5LCBwYXRoID0gaGVyZTo6aGVyZSgiZmlndXJlcyIpLCBmaWxlbmFtZSA9InBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX2NvbG9yYnl0cmVuZF9OWl9XQ1NJX1NFVVNfc3VtbWVyX05FVVNmYWxsb25seS5qcGciLCBoZWlnaHQgPSA2LCB3aWR0aCA9IDExLCB1bml0ID0gImluIikKCgoKYGBgCgpNZXJnZSBCQyB2ZXJzdXMgWWVhciBwbG90IHdpdGggR0FNUyBhbmQgUmVnaW9uIHZzLiBjb2VmZmljaWVudCBwbG90IGZvciBMTUVScwoKYGBge3J9CgpCQ190b3RhbF9HQU1fTE1FUl9tZXJnZV9sZWdlbmQgPC0gZ2dkcmF3KHhsaW0gPSBjKDAsIDQwLjUpLCB5bGltID0gYygwLCAyMSkpICsKICAgIGRyYXdfcGxvdChwb2ludHNfd2F2eWxpbmVzX2JyYXlfdG90YWxfeWVhcl9yZWR1Y2VkX2dhbSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB4ID0gMSwgeSA9IDEsIHdpZHRoID0gMjAsIGhlaWdodCA9IDIwKSArCiAgICBkcmF3X3Bsb3QoQkNfdG90YWxfRGlzc2ltaWxhcml0eV9Db2VmX2Vycm9yYmFyX3JlZHVjZWQgKwogICAgICAgIHRoZW1lKGxlZ2VuZC5rZXkuc2l6ZSA9IHVuaXQoMC41LCAnY20nKSwgI2NoYW5nZSBsZWdlbmQga2V5IHNpemUKICAgICAgICBsZWdlbmQudGl0bGUgPSBlbGVtZW50X3RleHQoc2l6ZT0xNiksICNjaGFuZ2UgbGVnZW5kIHRpdGxlIGZvbnQgc2l6ZQogICAgICAgIGxlZ2VuZC50ZXh0ID0gZWxlbWVudF90ZXh0KHNpemU9MTQpKSwgI2NoYW5nZSBsZWdlbmQgdGV4dCBmb250IHNpemUpLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHggPSAyMCwgeSA9IDEsIHdpZHRoID0xOSwgaGVpZ2h0ID0gMjApICsKICAgIGRyYXdfcGxvdChnZXRfbGVnZW5kKGRpcmVjdGlvbmFsX2NoYW5nZV9sZWdlbmRfcGxvdCArIAogICAgICB0aGVtZShsZWdlbmQua2V5LnNpemUgPSB1bml0KDAuNSwgJ2NtJyksICNjaGFuZ2UgbGVnZW5kIGtleSBzaXplCiAgICAgICAgbGVnZW5kLnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemU9MTUpLCAjY2hhbmdlIGxlZ2VuZCB0aXRsZSBmb250IHNpemUKICAgICAgICBsZWdlbmQudGV4dCA9IGVsZW1lbnRfdGV4dChzaXplPTEzKSkpLCAjY2hhbmdlIGxlZ2VuZCB0ZXh0IGZvbnQgc2l6ZSkKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgeCA9IDI3LCB5ID0gMTAsIHdpZHRoID0gMywgaGVpZ2h0ID0gMikgKwogIGdlb21fdGV4dChhZXMoeCA9IDIsIHkgPSAyMC43KSwgbGFiZWwgPSAoImEuIiksIHNpemUgPTgsIGZvbnRmYWNlID0gImJvbGQiKSArCiAgZ2VvbV90ZXh0KGFlcyh4ID0yMCwgeSA9IDIwLjcpLCBsYWJlbCA9ICgiYi4iKSwgc2l6ZSA9OCwgZm9udGZhY2UgPSAiYm9sZCIpCgoKZ2dzYXZlKEJDX3RvdGFsX0dBTV9MTUVSX21lcmdlX2xlZ2VuZCwgcGF0aCA9IGhlcmU6OmhlcmUoImZpZ3VyZXMiKSwgZmlsZW5hbWUgPSAiQkNfdG90YWxfR0FNX0xNRVJfbWVyZ2VfbGVnZW5kLnBuZyIsIGhlaWdodCA9IDgsIHdpZHRoID0gMTQsIHVuaXRzID0gImluIikKCiNBTFQgR1JFWSBTQ0FMRQpCQ190b3RhbF9HQU1fTE1FUl9tZXJnZV9sZWdlbmRfZ3JleXNjYWxlIDwtIGdnZHJhdyh4bGltID0gYygwLCA0MC41KSwgeWxpbSA9IGMoMCwgMjEpKSArCiAgICBkcmF3X3Bsb3QocG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1fZ3JleXNjYWxlLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHggPSAxLCB5ID0gMSwgd2lkdGggPSAyMCwgaGVpZ2h0ID0gMjApICsKICAgIGRyYXdfcGxvdChCQ190b3RhbF9EaXNzaW1pbGFyaXR5X0NvZWZfZXJyb3JiYXJfcmVkdWNlZF9ncmV5c2NhbGUgKwogICAgICAgIHRoZW1lKGxlZ2VuZC5rZXkuc2l6ZSA9IHVuaXQoMC41LCAnY20nKSwgI2NoYW5nZSBsZWdlbmQga2V5IHNpemUKICAgICAgICBsZWdlbmQudGl0bGUgPSBlbGVtZW50X3RleHQoc2l6ZT0xNiksICNjaGFuZ2UgbGVnZW5kIHRpdGxlIGZvbnQgc2l6ZQogICAgICAgIGxlZ2VuZC50ZXh0ID0gZWxlbWVudF90ZXh0KHNpemU9MTQpKSwgI2NoYW5nZSBsZWdlbmQgdGV4dCBmb250IHNpemUpLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHggPSAyMCwgeSA9IDEsIHdpZHRoID0gMTksIGhlaWdodCA9IDIwKSArCiAgICBkcmF3X3Bsb3QoZ2V0X2xlZ2VuZChkaXJlY3Rpb25hbF9jaGFuZ2VfbGVnZW5kX3Bsb3QgKyAKICAgICAgdGhlbWUobGVnZW5kLmtleS5zaXplID0gdW5pdCgwLjUsICdjbScpLCAjY2hhbmdlIGxlZ2VuZCBrZXkgc2l6ZQogICAgICAgIGxlZ2VuZC50aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplPTE1KSwgI2NoYW5nZSBsZWdlbmQgdGl0bGUgZm9udCBzaXplCiAgICAgICAgbGVnZW5kLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZT0xMykpKSwgI2NoYW5nZSBsZWdlbmQgdGV4dCBmb250IHNpemUpCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgeCA9IDI3LCB5ID0gMTAsIHdpZHRoID0gMywgaGVpZ2h0ID0gMikgKwogIGdlb21fdGV4dChhZXMoeCA9IDIsIHkgPSAyMC43KSwgbGFiZWwgPSAoImEuIiksIHNpemUgPTgsIGZvbnRmYWNlID0gImJvbGQiKSArCiAgZ2VvbV90ZXh0KGFlcyh4ID0yMCwgeSA9IDIwLjcpLCBsYWJlbCA9ICgiYi4iKSwgc2l6ZSA9OCwgZm9udGZhY2UgPSAiYm9sZCIpCgpnZ3NhdmUoQkNfdG90YWxfR0FNX0xNRVJfbWVyZ2VfbGVnZW5kX2dyZXlzY2FsZSwgcGF0aCA9IGhlcmU6OmhlcmUoImZpZ3VyZXMiKSwgZmlsZW5hbWUgPSAiQkNfdG90YWxfR0FNX0xNRVJfbWVyZ2VfbGVnZW5kX2dyZXlzY2FsZS5wbmciLCBoZWlnaHQgPSA4LCB3aWR0aCA9IDE0LCB1bml0cyA9ICJpbiIpCgojQUxUIENPTE9SIEJZIFRSRU5ECkJDX3RvdGFsX0dBTV9MTUVSX21lcmdlX2xlZ2VuZF9jb2xvcmJ5dHJlbmQgPC0gZ2dkcmF3KHhsaW0gPSBjKDAsIDQwLjUpLCB5bGltID0gYygwLCAyMSkpICsKICAgIGRyYXdfcGxvdChwb2ludHNfd2F2eWxpbmVzX2JyYXlfdG90YWxfeWVhcl9yZWR1Y2VkX2dhbV9jb2xvcmJ5dHJlbmQsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgeCA9IDEsIHkgPSAxLCB3aWR0aCA9IDIwLCBoZWlnaHQgPSAyMCkgKwogICAgZHJhd19wbG90KEJDX3RvdGFsX0Rpc3NpbWlsYXJpdHlfQ29lZl9lcnJvcmJhcl9yZWR1Y2VkX2NvbG9yYnl0cmVuZCArCiAgICAgICAgdGhlbWUobGVnZW5kLmtleS5zaXplID0gdW5pdCgwLjUsICdjbScpLCAjY2hhbmdlIGxlZ2VuZCBrZXkgc2l6ZQogICAgICAgIyBsZWdlbmQudGl0bGUgPSBlbGVtZW50X3RleHQoc2l6ZT0xNiksICNjaGFuZ2UgbGVnZW5kIHRpdGxlIGZvbnQgc2l6ZQogICAgICAgIyBsZWdlbmQudGV4dCA9IGVsZW1lbnRfdGV4dChzaXplPTE0KQogICAgICAgKSwgI2NoYW5nZSBsZWdlbmQgdGV4dCBmb250IHNpemUpLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHggPSAyMCwgeSA9IDEsIHdpZHRoID0gMTksIGhlaWdodCA9IDIwKSArCiAgICBkcmF3X3Bsb3QoZ2V0X2xlZ2VuZChkaXJlY3Rpb25hbF9jaGFuZ2VfbGVnZW5kX3Bsb3RfY29sb3JieXRyZW5kICsgCiAgICAgIHRoZW1lKGxlZ2VuZC5rZXkuc2l6ZSA9IHVuaXQoMC41LCAnY20nKSwgI2NoYW5nZSBsZWdlbmQga2V5IHNpemUKICAgICAgICBsZWdlbmQudGl0bGUgPSBlbGVtZW50X3RleHQoc2l6ZT0xNiksICNjaGFuZ2UgbGVnZW5kIHRpdGxlIGZvbnQgc2l6ZQogICAgICAgIGxlZ2VuZC50ZXh0ID0gZWxlbWVudF90ZXh0KHNpemU9MTUpKSksICNjaGFuZ2UgbGVnZW5kIHRleHQgZm9udCBzaXplKQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHggPSAyNywgeSA9IDgsIHdpZHRoID0gMywgaGVpZ2h0ID0gMikgKwogIGdlb21fdGV4dChhZXMoeCA9IDIsIHkgPSAyMC43KSwgbGFiZWwgPSAoImEuIiksIHNpemUgPTgsIGZvbnRmYWNlID0gImJvbGQiKSArCiAgZ2VvbV90ZXh0KGFlcyh4ID0yMCwgeSA9IDIwLjcpLCBsYWJlbCA9ICgiYi4iKSwgc2l6ZSA9OCwgZm9udGZhY2UgPSAiYm9sZCIpCgpnZ3NhdmUoQkNfdG90YWxfR0FNX0xNRVJfbWVyZ2VfbGVnZW5kX2NvbG9yYnl0cmVuZCwgcGF0aCA9IGhlcmU6OmhlcmUoImZpZ3VyZXMiKSwgZmlsZW5hbWUgPSAiQkNfdG90YWxfR0FNX0xNRVJfbWVyZ2VfbGVnZW5kX2NvbG9yYnl0cmVuZC5wbmciLCBoZWlnaHQgPSA4LCB3aWR0aCA9IDE0LCB1bml0cyA9ICJpbiIpCgpgYGAKCklzIGVhY2ggc3VydmV5IGJldHRlciBkZXNjcmliZWQgYnkgbG1lciBvciBieSBHQU0/CmBgYHtyfQoKIyBGdW5jdGlvbiB0byBjb21wYXJlIEdBTSBhbmQgbGluZWFyIG1vZGVsIGZpdHMKY29tcGFyZV9tb2RlbHMgPC0gZnVuY3Rpb24oc3Vic2V0X2RhdGEpIHsKICAjIEZpdCBHQU0gbW9kZWwgd2l0aCAiZnMiIGJhc2lzIGFuZCBzaW5nbGUgc21vb3RoIHRlcm0KICBnYW1fbW9kZWwgPC0gZ2FtKGJyYXlfY3VydGlzX2Rpc3NpbWlsYXJpdHlfdG90YWxfbWVhbiB+IHMoeWVhciwgYnMgPSAiZnMiLCBtID0gMSksIGRhdGEgPSBzdWJzZXRfZGF0YSkKICAKICAjIEZpdCBsaW5lYXIgbW9kZWwKICBsbV9tb2RlbCA8LSBsbShicmF5X2N1cnRpc19kaXNzaW1pbGFyaXR5X3RvdGFsX21lYW4gfiB5ZWFyLCBkYXRhID0gc3Vic2V0X2RhdGEpCiAgCiAgIyBDYWxjdWxhdGUgQUlDYyB2YWx1ZXMKICBhaWNjX2dhbSA8LSBBSUNjKGdhbV9tb2RlbCkKICBhaWNjX2xtIDwtIEFJQ2MobG1fbW9kZWwpCiAgCiAgY29tcGFyaXNvbiA8LSBpZmVsc2UoKGFicyhhaWNjX2dhbS1haWNjX2xtKSA8IDIgfCAoYWljY19nYW0gPT0gYWljY19sbSkpLCJTYW1lIiwKICAgICAgICAgICAgICAgICAgICAgICBpZmVsc2UoKChhYnMoYWljY19nYW0tYWljY19sbSkgPiAyICYgKGFpY2NfZ2FtIDwgYWljY19sbSkpKSwiR0FNIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgaWZlbHNlKCgoYWJzKGFpY2NfZ2FtLWFpY2NfbG0pID4gMiAmIChhaWNjX2dhbSA+IGFpY2NfbG0pKSksICJMaW5lYXIiLE5BKSkpCiAgCiAgIHJldHVybihkYXRhLnRhYmxlKFN1cnZleV9OYW1lX1NlYXNvbiA9IHVuaXF1ZShzdWJzZXRfZGF0YSRTdXJ2ZXlfTmFtZV9TZWFzb24pWzFdLAogICAgICAgICAgICAgICAgICAgIGNvbXBhcmlzb24gPSBjb21wYXJpc29uLAogICAgICAgICAgICAgICAgICAgIGFpY2NfZ2FtID0gYWljY19nYW0sCiAgICAgICAgICAgICAgICAgICAgYWljY19sbSA9IGFpY2NfbG0pKQp9CgpHQU1fTE1fbW9kZWxfY29tcGFyaXNvbnMgPC0gZGF0YS50YWJsZSgpCgojIENvbXBhcmUgbW9kZWxzIGZvciBlYWNoIHN1cnZleSBuYW1lIGFuZCBiaW5kIHRoZSByZXN1bHRzCmZvciAoaSBpbiAxOmxlbmd0aChhbGxfc3VydmV5X25hbWVzKSkgewogIEdBTV9MTV9tb2RlbF9jb21wYXJpc29uc19zaW5nbGUgPC0gY29tcGFyZV9tb2RlbHMoZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yW1N1cnZleV9OYW1lX1NlYXNvbiA9PSBhbGxfc3VydmV5X25hbWVzW2ldXSkKICBHQU1fTE1fbW9kZWxfY29tcGFyaXNvbnMgPC0gcmJpbmQoR0FNX0xNX21vZGVsX2NvbXBhcmlzb25zLCBHQU1fTE1fbW9kZWxfY29tcGFyaXNvbnNfc2luZ2xlKQp9CgpHQU1fTE1fbW9kZWxfY29tcGFyaXNvbnMKCiNyb3VuZCB0byAyIHNpZ25pZmljYW50IGZpZ3VyZXMKCiNuYW1lcyBvZiBudW1lcmljIGNvbHVtbnMKbnVtZXJpY19jb2xzIDwtIG5hbWVzKEdBTV9MTV9tb2RlbF9jb21wYXJpc29ucylbc2FwcGx5KEdBTV9MTV9tb2RlbF9jb21wYXJpc29ucywgaXMubnVtZXJpYyldCgojIEFwcGx5IHNpZ25pZigpIG9ubHkgdG8gbnVtZXJpYyBjb2x1bW5zCkdBTV9MTV9tb2RlbF9jb21wYXJpc29uc1ssIChudW1lcmljX2NvbHMpIDo9IGxhcHBseSguU0QsIGZ1bmN0aW9uKHgpIGlmIChpcy5udW1lcmljKHgpKSBzaWduaWYoeCwgZGlnaXRzID0gMykgZWxzZSB4KSwgLlNEY29scyA9IG51bWVyaWNfY29sc10KCgpWaWV3KEdBTV9MTV9tb2RlbF9jb21wYXJpc29ucykKCmNvbG5hbWVzKEdBTV9MTV9tb2RlbF9jb21wYXJpc29ucykgPC0gYygiU3VydmV5IiwiTW9kZWwgY29tcGFyaXNvbiIsICJHQU0gQUlDYyIsIkxNIEFJQ2MiKQoKZndyaXRlKEdBTV9MTV9tb2RlbF9jb21wYXJpc29ucywgaGVyZTo6aGVyZSgib3V0cHV0IiwiR0FNX0xNX21vZGVsX2NvbXBhcmlzb25zLmNzdiIpKQpgYGAKCg==